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Abstract 


Statistical tools, previously developed for nonlinear least-squares 
estimation of multivariate sensor calibration parameters and the 
associated calibration uncertainty analysis, have been applied to 
single- and multiple-axis inertial model attitude sensors used in wind 
tunnel testing to measure angle of attack and roll angle. The anal- 
ysis provides confidence and prediction intervals of calibrated sensor 
measurement uncertainty as functions of applied input, pitch and roll 
angles. A comparative performance study of various experimental 
designs for inertial sensor calibration is presented along with corrob- 
orating experimental data. The importance of replicated calibrations 
over extended time periods has been emphasized; replication provides 
independent estimates of calibration precision and bias uncertainties, 
statistical tests for calibration or modeling bias uncertainty, and sta- 
tistical tests for sensor parameter drift over time. A set of recom- 
mendations for a new standardized model attitude sensor calibration 
method and usage procedures is included. The statistical information 
provided by these procedures is necessary for the uncertainty analy- 
sis of aerospace test results now required by industrial users of wind 
tunnel test facilities. 

1. Introduction 

The standard instrumentation used at the Langley Research Center (LaRC) for measuring 
model attitude in the wind tunnel is the inertial angle of attack (AOA) sensor package described 
in reference 1. Langley Research Center has employed the inertial sensor as the primary AOA 
measurement system during the past 30 years. Various aspects of inertial model attitude 
measurement have been subsequently reported in references 2 to 4. In particular, reference 2 
describes data reduction techniques for model attitude measurements in pitch and roll and pitch 
measurement only at zero roll. Typically, the LaRC AOA package provides static model attitude 
measurements at accuracies of ±0.01°. 

Because of signal-to-noise ratios as low as — 100 dB commonly encountered in wind tunnel test 
facilities, heavy low-pass filtering in the bandwidth range of 0.3 to 0.6 Hz is necessary for static 
attitude measurement (ref. 3). Therefore the inertial system is suitable only as a static attitude 
measurement device and is not useful for dynamic attitude measurement. In addition, the inertial 
accelerometer has been found to exhibit an offset error due to centrifugal forces developed in the 
presence of repetitive model motion in yaw and pitch encountered at high dynamic levels during 
tests, as discussed in reference 4. Although optical sensors, which are insensitive to centrifugal 
errors, are used increasingly for both static and dynamic model attitude measurement, the 
inertial sensor remains important for high-precision primary measurement, calibration of optical 
systems, and optical system backup during poor test section visibility. 

Inertial model attitude sensor packages have been calibrated at LaRC by means of four- and 
six-point tumble tests. The tumble test technique, easy to implement through the use of simple 
precision leveling devices, has been adequate in the past. It, however, does not provide adequate 
spatial resolution for modeling precision or statistical uncertainty information now required by 
test facility users. Also, current calibration procedures do not employ replication, necessary for 
independent estimation of sensor bias and precision uncertainties and for assessment of long-term 
drift. 

Multiple-point replicated calibration is now feasible and convenient through use of the 
automatically controlled calibration dividing head and modern computerized control and data 



acquisition systems. Statistical tools recently developed in reference 5 for general estimation of 
multivariate sensor calibration parameters and the associated calibration uncertainty analysis 
are applied in this publication to multiple- point replicated calibration of inertial AOA packages. 
These statistical tools, applied to one-, two-, and three-axis inertial sensor packages, allow 
comparison of experimental designs for calibration, computation of calibration confidence 
intervals, and prediction intervals as functions of applied inputs, independent estimation 
of calibration bias and precision uncertainties, and detection of long-term parameter drift. 
Experimental calibration data are presented to demonstrate and verify the efficacy of the 
technique . 

Based on the theoretical analysis and experimental calibration results, a set of recommenda- 
tions for model attitude sensor calibration and usage is proposed. The recommended procedures 
may be readily implemented by means of modern automated calibration apparatus. The sta- 
tistical information thus provided, not previously available to test facility users, is necessary 
for determination of overall uncertainty of aerospace test results now required by industrial test 
facility users. 

2. Pitch and Roll Measurement 

2.1. Angle of Attack Measurement Without Roll 

Use of the single-axis inertial angle of attack (AOA) sensor in wind tunnel facilities without 
roll allows simplified data reduction, as described in reference 2; the uncertainty analysis 
described briefly in reference 5 is extended here. Misalignment of the accelerometer sensitive 
axis with respect to the AOA package x-axis is represented by the angle, denoted by <j> a , between 
the projection of the sensitive axis onto the x-z (pitch) plane and the x-axis. Roll angles during 
calibration and facility usage are assumed to remain zero. The sensor output is given by the 
following equation: 

v a = b a + S a sin (a - <j > a ) (1) 

where v a is the sensor output in volts, b a is the sensor offset in volts, S a is the sensitivity in volts 
per g unit, a is the pitch angle in radians, and <j> a is the pitch misalignment angle in radians. 
Note that acceleration of gravity g is normalized to unity in all equations. 

2.2. Model Attitude Measurement With Roll 

For single-axis or multiple-axis attitude measurement with roll, the inertial sensor axis 
misalignment must be characterized in three-dimensional (3-D) space. At LaRC the sensitive 
axis of the x-axis sensor is represented as lying on the surface of a cone, aligned with the x-axis 
of the sensor package, whose vertex is located at the origin of the package coordinate system. 
The semivertex angle of the cone, denoted by Q, is termed the “coning angle.” Looking in the 
positive x direction, the angular position of the pitch sensor axis on the surface of the cone is 
specified by angle A x , measured counterclockwise from the positive y-axis to the pitch sensor 
axis; angle A x is termed the “azimuth angle.” As indicated in appendix A and reference 2, the 
sensor output equation is given by the following form: 


v x = b x + Sjjcos fU sin a — sin Q x cos a sin(R + Aj.)] 


( 2 ) 


where R denotes roll angle and subscript x denotes pitch sensor parameters. Angles are in 
radians. If roll angle R is known, input angle a is inferred by inverting equation (2) to obtain 


a = arcs in 


( 'Ux bx) j Sx 

1 ^/cos 2 Q x + sin' 1 (R + A,,) siiU Q x 


+ arctan [tan Q x sin (R + A x )\ 


(3) 
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Multiple-axis inertial attitude measurement packages, designed for simultaneous measure- 
ment of pitch and roll angles, employ two orthogonally placed accelerometers aligned nomi- 
nally with the x- and y- axes of the model, or three orthogonally placed accelerometers aligned 
nominally with the x -, y-, and z-axes of the model. Coning angles fi v and Q, and azimuth angles 
A y and A z for the y-axis and z-axis sensors are defined analogously to Q x and A x . The s-axis 
sensor output is given by equation (2). The y-axis sensor output, obtained in appendix A, is 
found to be 

v y — by — S' v [cos sin R cos a — sin Q y (sin A,, sin a — cos A y cos R cos a)] (4) 

Given observed outputs v x and v y , the corresponding inputs a and R are inferred by 

simultaneous solution of equations (2) and (4) via an iterative method. However, as shown 
later a useful solution does not exist near a = ±90° or R = ±90°, where the 2x2 Jacobian 
matrix of the system of equations (2) and (4) with respect to a and R becomes singular or poorly 
conditioned. It can be shown that the Jacobian matrix must be nonsingular for the existence of 
a solution (ref. 6). 

As shown later, the singularities near R = ±90° are eliminated by addition of the z-axis 
sensor, whose output, obtained in appendix A, is found to be 

v z = b z — S z [cos cos R cos a — sin Q, (cos A z sin a — sin A z sin R cos a)] (5) 

The 3 x 2 J acobian matrix of the system of equations (2), (4), and (5) has rank 1 at a = ±90°, 

and rank 2 elsewhere for f2 < 10° as is shown subsequently. Inputs a and R are estimated by 
least-squares solution of the overdetermined system of equations (2), (4), and (5), provided that 
the Jacobian matrix has rank 2. At a = ±90°, estimated pitch angle can be determined within 
the accuracy of the y-axis and z-axis sensors, although roll angle cannot be determined. Note 
that calibration parameters 6, S, f2, and A of sensors x, y, and z are independently determined. 

3. Experimental Designs for Calibration 

Experimental designs for calibration of the single-axis AOA sensor without roll, the single- 
axis pitch sensor with roll, and the multiple-axis package are now analyzed by using nonlinear 
multivariate uncertainty analysis techniques and notation developed in reference 5. Let c a denote 
the 3x1 parameter vector for the single-axis sensor without roll as follows: 

Ca = [be So (paY (6) 

and let z denote the vector of independent variables, which contains the single element a. The 
calibration experimental design D consists of /f-elertienl, set fi Q , = {aq,...,a A -} C [a m i n , a mox ]. 
The K x 1 design matrix Z is then 

Z = [a 1 ...a K ] r (7) 

Similarly, let c x , c y , and c z denote 4x1 vectors of x-, y-, and z-axis sensor parameters with roll; 
therefore, 

c* = [b x S x Q x A x f ' 
c y = [ b u S y O v AyY > 
c = [b S z L> AYY , 

and let z denote the 1x2 vector of independent variables 

z = [a R] (9) 


( 8 ) 
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The calibration experimental design contains K pitch-roll angle pairs, where the pitch angle is 
selected from set fi Q C [« mnn a max ] containing N values, and the roll angle is selected from set 
B^. Q [-Rmim -RmJ containing M values; thus 


6 n = {a i ,a - 2 , , a ,v} 

Bfl = {Ah, R-2, ■ ■ ■ , Rm} 


( 10 ) 


The experimental design of primary interest, denoted by D, contains K = MN ordered pairs 
from sets fi Q and Bj?, represented by K x 2 design matrix Z as 


Oil 

a r ■ 

■ ■ a L 

“2 

Oi-2 ■ 

• • a 2 

a , v 


Ri 

R-2 ■ 

Rm 

Ry 

R-2 • 

' ' Rm ‘ ■ 

• Ry 

R-2 


Qy 

R> 


M 


(11) 


Although, as is shown, design D has desirable properties, its possibly large cardinality may 
become experimentally impractical. Fractional experimental designs constructed as subsets of D 
are described later and provide more efficient calibration with adequate prediction uncertainties. 
The considerable available literature on design of efficient experiments is not reviewed in this 
publication. 

Let the corresponding x-, y-, and z-axis sensor output observations be denoted by K x 1 
vectors v„ v y , and v z as follows: 


V* = 

[V X 1 

V x -2 ■ 

■ ■ v xK 

f' 

v .v = 

[ Vyl 

Vy2 ■ ■ 

■ VyK 

I" ^ 

V ■/ = 

[ Vzl 

V :-2 ■ ■ 

■ VzK. 

w 


(12) 


3.1. Observed Sensor Outputs 

At the k th calibration point of the single-axis sensor without roll, where k = 1 , ... , K , 
element k of observation vector v Q is obtained from equation (1) as 


v ak = f a {c a ,z k ) = b Q + 5), [cos <j> a sin a k - sin <j> a cos a*] (13) 

Similarly, for sensors with roll, vectors v t , v v , and \ z are obtained by using equations (2), (4), 
and (5) as 


v Xl: = /*(c*,z*) = b x + S x [cos Q a . sin a k - sin cos a k sin (R k + A,)] (14) 

= by — S y [cos Q„ sin R k cos a k — sin O v (sin A y sin a k — cos A y cos R k cos a^)] (15) 

and 


*L { . = fz(Cz,Zk) 

= b, — S ,. [cos S2, cos R k cos a k — sin Q, (cos A z sin a k — sin A z sin R k cos a 4 )] (16) 


Note that equation (13) is a special case of equation (14), where 


= <j> a ' 
A x = 7t/2 > 

Rk = 0 


(17) 
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For K observations, equations (13) to (16) are extended to vector function notation as 


^ x ^(^£>21) \f x (c x , Zj ) f x ic x , Z 2 ) ••• /.(c^Za-)] 


Vectors v Q , v„, and v z are defined analogously. 


3.2. Evaluation of Gradient Matrices 


The 3x1 gradient vector of f a (c a ,z) with respect to c a is given by 


f _ df Q (c Q ,z) _ r d/ n (c a ,z) df a (c a ,z) d/ Q (c a ,z) 

° c — dc a ~ db a dS a d<j> a 


The 4x1 gradient vectors — /„( c„,z) with respect to c x , f y (c y ,z) with respect to c„, and c z ,z) 
with respect to c z — are obtained as follows: 


df(c,z) _ [ df(c,z) df( c,z) df{ c,z) df{c, z) 
dc ~ db dS dQ dA 


Element-by-element evaluation of equation (20) for the a: -axis sensor is as follows: 

_ df x { c„z) _ 

1x1 _ db x 

f df x (c x ,z) O' o ! 1 A \ 

j x 9 = = cos sin a — sin \l e cos a sin (it + A x ) 

dS x 

fxQ = ^ = ->5* [sin Q.v sin a + cos O,. cos a sin (R + A x )] = S x (f>xQ 

fxA = = -S x sin Cl x cos a cos (R + A x ) = S x w x <j> xA 


where 


<j> x Q = —sin sin a — cos cos a sin ( R + A x ) 
<j> x _4 = —cos a cos (R+ A x ) 
w r = sin 


To evaluate the gradient terms of equation (19) for the single-axis sensor without roll, substitute 
the values of equation (17) into equations (21) to (23). 


Similarly equation (20) is evaluated for the y- 


axis sensor 


as follows: 


dfJ c„,z) 

u = ; )h = 1 ( 28 ) 

dfy(c y ,z) 

fyc. = — = —cos O y cos a sin R + sin Q v (sin A,, sin a — cos A y cos a cos R ) (29) 

ob„ 


= df u { Cy,z) 

ha dQ„ 


= S'jjsin Qy cos a sin R + cos fi y (sin A y sin a — cos A y cos a cos R)] 


= Sy(f>y(l 
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_ dfy(Cy,Z) 


f Vi = —jj ^ — = S y sin Q y (cos A y sin a + sin A y cos R cos a) = S y w y <f> yA (31) 


<j>y q = sin Qy cos a sin R + cos Q v (sin A y sin a — cos A y cos a cos R) 
(j > y 4 = cos Ay sin a + sin A y cos R cos a 


w y — sin Q y 


Evaluation of equation (20) is similar for the z-axis sensor as follows: 

/„ S = 1 (35) 

db z 

f z 9 = — - = —cos Q z cos a cos R + sin Q ,(cos A z sin a — sin A z cos a sin R) (36) 
dS z 

f z Q = ’ - = 5 ? [sin Qy cos a cos R + cos Q;.(cos A z sin a — sin A z cos a sin R)] 

= S,K (37) 

f z 4 = 4 ’ - — —S z sin Q z (sin A., sin a + cos A z sin R cos a) = S z w z <j > z 4 (38) 

oA z 

where 

<j) z o. = sin Qj cos a cos i? + cos Q ; (cos A sin a — sin A z cos a sin R) (39) 

(f> zA = — (sin A z sin a + cos A z sin R cos a) (40) 

w z = sin Q z (41) 

For calibration of sensor packages with roll, define K x 4 gradient matrices F lc , F Vc , and F Zc , 
obtained from equation ( 20 ) as 


df(c,Z) . . . I 

Fc — ^ -= = ; f 9 ! f q ; fA 


where f f?, f^, and f 4 denote columns 1, 2, 3, and 4, respectively, of matrix F c . The K x 3 
matrix F QC is similarly defined for the single-axis sensor without roll. 

Reference 5 shows that the leashsquares estimate of c, denoted by c, is individually obtained 
for sensor x, y, or z by solving the following K x 1 system of nonlinear equations for c: 


ITtt-1 I M , 


h(v,c)= [v-f(c,Z)] i U Y 1 — (c,Z) =0 


where v is the K x 1 vector of observed outputs, and ct 2 e \Jy is the K x K output uncertainty 
covariance matrix, where <r| is the measurement variance. The Lx L moment matrix R (ref. 5) 
is given by the following equation: 


RE ^ = F ' UYlFt + H£ 
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where the ij th element of L x L matrix H £ is given by 


h Hj = [v -f(c,Z)] T U Y 1 f cc .. (45) 

where f CC; . is the ijth column vector of length K contained in L x L x K array F cc defined by 

where 1 < i, j < L, L = 3 without roll; and L = 4 with roll. Matrix H^, evaluated in 
appendix B, is negligible unless the least-squares residuals are large. Indeed, note that term 
e = [v — f(c, Z)] in equation (45) equals the vector of residuals following least-squares estimation 
of c. The norm of e, equal to the root sum of squares of its elements, is defined as 


' A’ \ ^ 

=(e T e)W (47) 

,k=l / 

Reference 5 shows that the expected value of ||e|| equals (K— L)F 2 a E , where a E is the standard 
deviation of the measurement error. Therefore, if a E is small, matrix can be neglected in 
equation (44) for uncertainty analysis. See appendix B for details. 

The standard error Se, defined individually for sensor x , y , or _ as 

c _ (e T e ) 1/2 

p — ■ = 

s/I< - L 

provides an unbiased estimate of a E . For the special case where U y = I and where H s can be 
neglected, moment matrix R becomes 




R = F^F c (49) 

The covariance matrix of estimated parameter vector c is then given by (ref. 5) 

£ c = 4R- 1 (50) 


A confidence ellipsoid for c at confidence level 1 —cr is defined by the following inequality (ref. 5): 

(c - c) t R~ j (c - c) < ( K - L)S 2 E F LK _ L (a) (51) 

where F L K _ L ( a) is the a- percentile value of the F-distribution with L, K — L degrees of freedom. 


3.3. Sensor Output Variance Function 

In reference 5, the variance function cr 2 (z) of predicted outputs v x , v u , and v zi respectively, 
for sensor x, y , and z is given by the following quadratic form: 


oj(z) _ f r T (z)S c f c (z) 




f c T (z)R %(«) 


(52) 


The following three theorems, proved in appendix C, show that the output variance functions 
of the x-, y-, and z-axis sensors are independent of the corresponding parameter vector c for 
any calibration experimental design. 
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Theorem I: Sensor output variance function <x 2 (z) is independent of calibration parameters b 
and S 

Theorem II: Sensor output variance function <x 2 (z) is independent of calibration parameter Q 
Theorem III: S ensor output variance function cr(z) is independent of calibration parameter A 

Note in equation (127), proof of Theorem 1 in appendix C, that variance function cb 2 (z) is 
well-defined whenever matrix P (eq. (129)) is nonsingular. Thus <r 2 ( z) exists for w = sin 0 = 0 
where matrix R is singular. Matrix R is evaluated analytically in appendix D. 

From Theorems I to III, the conclusion is drawn that variance function <r’ 2 ( z) of predicted 
output v is independent of calibration parameters 6, S, Q, and A for the x-, y and 2 -axis 
sensors. Hence, sensor output uncertainty depends only upon experimental design values of a 
and R and measurement variance cr 2 E . 

3.4. Experimental Design Figure of Merit 

Box (ref. 7) defines a figure of merit V for any experimental design as the mean value of the 
output variance function over test volume Sr, normalized by the number of calibration points 
and the measurement variance. (See also ref. 5.) The value of V for experimental design D is 
obtained with the help of equation (147) as 


I< 


V = 


Pftmax rR max 

MN / / J.TD-1J.T. 

■'ttillir. J Ria 


:(z ) dx MN 


4> C P l 4 dR da 


/ dx 

r ft max fi 

/a 

d ft min 


(53) 


dR da 


Design figures of merit are equal for x-, y - , and a- ax is sensor output uncertainties. The 
numerator of equation (53), which contains integrals of cross products of the elements of gradient 
vector cf> c , is evaluated in appendix E as 


V N = / q R ( z) dx 


r Q max r i 
J cv.-.-.;.! J 


4> C V dRda = 

i —1 j — 1 


(54) 


where PD is the fjth element of the inverse of matrix P defined in equation (129) and terms 
1^.. are defined in appendix E. The figure-of-merit expression 

MNV n 
hi 

is obtained in appendix E. Definite integral I hi is defined for the £-axis sensor in equation (214). 
Values of V for selected experimental designs are given later. 

4. Evaluation of Variance Function for Special Experimental Designs 

4.1. Experimental Designs 

Three special calibration experimental designs, denoted by Do, Di, and T, are considered as 
follows: 

Minimal design D 0 : A special case of design D 

1. Pitch angle set fi Q contains N points in the closed interval [a lrim , a max ] 

2. Roll angle set 6^. contains M unique principal angle valued points, uniformly distributed 
over closed interval [—7 r, 7r — AT?], where A R = 2tt/M 
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Minimal design Dp A special case of design D 0 

1. Pitch angle set contains N unique principal angle valued points uniformly distributed 
and centered about zero over the closed interval [— a max , cr max ], although a max may equal 
7 r, where Act — 2a max /(AT — 1). 

2. Roll angle set equals that of design D 0 

Parts 1 of designs D 0 and Dj apply for calibration without roll. Designs D 0 and Dj may also 
be constructed of multiple copies of a minimal Do or Di design, respectively. For example, a 
typical pitch calibration proceeds from to followed the same points in reverse order 

from a max to The properties of design D variance functions derived in sections 4.2 and 4.3 
are preserved under reordering, randomization, and replication. 

Design T 

1. Six-point “tumble” calibration with roll 

The single-axis or multiple-axis sensor package with roll is calibrated only at cardinal 
angles; experimental design matrix Z is as follows: 


Z = 



0 0 


7T 

2 

0 


TT 

0 


0 0 ■ 

_7T 7T 

'2 2 _ 


(56) 


2. Four-point tumble calibration without roll 

The single-axis sensor package without roll is calibrated only at cardinal angles; experi- 
mental design matrix Z is as follows: 


Z = 



(57) 


Moment matrix R and its related matrix P are evaluated analytically in appendix D in 
equations (206) to (2 13) for computation of variance function <r 2 (z). Because cr^z) is independent 
of parameters b , S, Q, and A , the following parameter values are chosen for simplification: 


b = n = A = 

S = 1 



(58) 


The values listed in equations (17) are selected for computation of variance function <r 2 (z) 
without roll. 


4.2. Variance Function for Design D 0 

Sensor output variance cr^z*,) for design D 0 depends only on the number of pitch calibration 
points N, the number of roll calibration points M , the pitch angle calibration range a max , and 
the pitch angle a k as shown by the following. The output variance for x-, y-, and z-axis sensors 
is given by equation (147) as 

^*) = ^(z t )p-V e (z*) (59) 

a E 

where <f> c is defined in equation (128) and matrix P is evaluated in appendix D (eq. (129)). The 
following theorem, proved in appendix C, shows, for calibration with roll, that the pitch angle 
sensor output uncertainty is independent of roll angle R for design Do. 
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Theorem IV: Let roll angle calibration set B fi , defined in equation (10), contain K = NM points 
uniformly spaced over the interval [ — 7T , 7 r— A R], where M and N are integers, 
A R = and the principal value of each angle contained in fi^ occurs with 

the same frequency; then the pitch sensor output variance is independent of roll 
angle R. 

For calibration without roll, equations (21), (22), and (25), evaluated by using the parameter 
values of equations (17), become 

Lb = i 


fxS = sin a 

<f>xn = -c°s a . 


With the help of equations (176) to (203), 


(60) 



V 

,Q bb 

W, 

Pa bSl 

p, = 

V 

1 Q bS 

r «9.9 

Pa. 9f2 


Pa bQ 

Paso 

Pan n 

r «bb 

= N 


1 


(61) 


(62) 


where 


ra bS = Sa 

P<*bQ = — 

r ass = \{N-C 2a ) 

__1 Q 

P °. sn 2'- >2a 

Pa QQ = 2 ^ iV + C ' 2a ' > 
and where Sa, C a , S 2 a , an d C 2 a are defined in equations (165) and (166). 

4.3. Variance Function for Design Dj 

For design D^, matrix P for the x-, y-, and z-axis sensors simplifies to the following diagonal 
form for calibration with roll: 


rn 0 
0 r ss 
0 
0 


0 
0 

0 Pm 

0 0 


0 

0 

0 

pAA . 


(63) 


Inverse matrix P 1 is given by 


1 

r bb 


-L- o 


i 

rss 

0 

0 


0 

0 

1 

Pm 

0 


0 
0 
0 

J_ 

Paa J 


(64) 
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Combine equations (59) and (64) to obtain x-, y-, and z-axis sensor output variances as 


^O) _ 1 

•i I 


t | <t>A 


Equation (65) is evaluated for design D] with the help of equations (148) and equations (176) 
to (203); after simplification the normalized x-axis sensor variance is obtained as 


2 [At + C 2 a + (AT — 3C 2a ) cos 2 a] 
M(AC - Cl) 


where C 2q is defined in appendix D (eqs. (172)). It is shown in appendix D that 
cr (z) = cr vz {z) = cr vx ( z). Equation (66) shows that the variation of a 2 l:x (z k ) with aq, is concave 
upward about zero pitch for C- 2 a > AT/3 and concave downward about zero pitch for C 2a < AT/ 3 . 
Normally, maximum attitude measurement accuracy is desired near zero pitch. 

For calibration without roll via design Dj, variables S, A = 0 and S 2 a = 0; equations (62) 
change accordingly. The variance function is shown to be given by 


<r^(z) (l/2)(A r + C 2a ) — '2C a cos a + N cos 2 a 2 sin 2 a 

' 1 1 / — ) a ! VtC„)-r:; + at - c 2a 


(67) 


4.4. Variance Function for Design T 

For single- axis or multiple-axis six-point tumble calibration with roll, matrix P (eq. (129)) 
simplifies to the following diagonal form for x-, y - , and z-axis sensors: 

6 0 0 0 
0 2 0 0 
0 0 2 0 
0 0 0 2 

From equations (154) and (65), the variance function is 

• 7 ;’ It-' _ 2 

°i 3 

After multiplying by the number of calibration points, the normalized standard deviation is 
found to be equal to 2. 

For single-axis four-point tumble calibration without roll, matrix P becomes 

'4 0 O' 

P = 0 2 0 (70) 

.0 0 2 _ 

The variance function is 

°j a ( z ) = 3 
4 

After multiplying the variance function by the number of calibration points, the normalized 
standard deviation is found to be equal to 3 1//2 . 
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5. Confidence and Prediction Intervals 

5.1. Multiple- Axis Sensor Uncertainty 

For arbitrary input zo, the calibration confidence interval of the corresponding predicted 
sensor output v (h for sensor x, y , or z, is defined by the following expression: 

6v 0 = |uo - ro| < tir-4(a)5 B [fJ(zo)R^ 1 fc(zo)] 1/2 (72) 

where Se is the standard error of the regression and tj.(a) is the a-percentile value of the 
two-tailed t-distribution with k degrees of freedom, denoted the precision index (ref. 8); x, 
y , and z subscripts are elided. The corresponding prediction interval (ref. 5) of a single new 
measured output is defined as 


& V P B < Cc-4(“)S£ 


f c T (z 0 )R- 1 f c (z n ) + ^ 


1/2 


(73) 


where <r 2 is the variance of the new measurement and is the calibration measurement variance. 


5.2. Single-Axis Pitch Sensor Uncertainty With Roll 

New measurement data reduction for the single-axis pitch sensor with roll requires inde- 
pendent measurement of roll angle R whose variance, denoted by cr 2 R , is independent of the 
calibration uncertainties and the pitch sensor output measurement uncertainty. The calibration 
confidence interval is given by equation (72). The prediction interval is given by 


SvPti 5: f K-i{ a )SE 


fJ(z 0 )R- 1 f c (z 0 ) + — ~ — 


I 1/2 


where 

fxR - = S x sin Q x cos a n cos (R a + A x ) 

and where a a and R a are the new pitch and roll angles, respectively. 


(74) 


(75) 


5.3. Parametric Studies of Experimental Designs 

Figures 1 to 5 illustrate the variation of sensor output uncertainty with pitch at selected 
parameter values for various experimental designs. Recall that uncertainties for x-, y-, and z-axis 
sensor output are identical. Uncertainties are shown as standard deviation functions normalized 
by sensor measurement uncertainty and (M7V) 1 / 2 , where M and N are the number of roll 
and pitch calibration points, respectively. Note that calibrations without roll are normalized 
by N 1 !' 2 . Confidence intervals are readily obtained from normalized standard deviation curves. 
For comparison, normalized tumble test uncertainty curves are shown with those of the higher 
order experimental design in each of figures 1 and 2. Note that the low cardinality of 
tumble calibrations causes high calibration uncertainties compared with higher order calibration. 
Although the normalized tumble calibration uncertainties are comparable with those of the 
higher order designs, the unnormalized tumble calibration uncertainties will increase by the 
factor (65/4) 1 / 2 in figure 1 and by (65/6) 1 / 2 in figure 2 compared with the uncertainties of the 
higher order designs. 

For comparison, table 1 presents the normalized mean standard deviations V 1 !' 2 , where V 
is the figure of merit defined in equation (53), for calibration designs with roll from figures 2 
to 5, evaluated over the calibration range. In addition, the normalized mean standard deviations 
evaluated over reduced usage ranges, denoted by V R ' 2 , are shown. 
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5. 3. 1 Single-axis pitch sensor without roll. Figure 1 illustrates the variation of normalized 
sensor output standard deviation with pitch angle for design Dj for calibrations over the ranges 
from —30° to 30°, —45° to 45°, —90° to 90°, and —180° to 180°, respectively, for N = 65. The 
constant normalized standard deviation for the four-point tumble calibration is shown in each 
figure for comparison. Note in figures 1(a), (b), and (c) that sensor uncertainty is low within 
the center 50 percent of the calibration range and increases rapidly outside the center range. 
Calibration from —180° to 180° produces nearly constant uncertainty approximately equal to 
that for the four-point tumble calibration and at a level 17 percent greater than that in the 
center ranges of the calibration designs from —90° to 90° and less. 

5.3.2. Single- or multiple-axis attitude sensor with roll. Some effects of spacing test 
points uniformly and nonuniformly on the mean normalized standard deviation using designs D 0 
and Dj are illustrated in figures 2 through 5 and summarized in table 1. Figure 2 illustrates the 
variation of sensor output standard deviation a vx with pitch angle for design Dj , for maximum 
pitch calibration angles of 30°, 45°, 90°, and 180°, respectively, and for values of N from the set 
{5, 9, 17, 33, 65}. For comparison, the constant normalized standard deviation for the six-point 
tumble calibration design T is indicated in each figure. As shown in Theorem IV in section 4.2, 
cr vx is independent of roll with design D 0 and, hence, with design Dj. From equation (66), the 
normalized curves are independent of M . Note in figures 2(a) and (b) that the uncertainty 
curves concave upward about 0° for calibration designs with a max < 45°. 

Figures 2(c) and (d) show that calibration for a = —90° to 90° and —180° to 180° produce 
uncertainty curves concaved downward about, 0° with significantly greater uncertainty at 0° than 
at ±90°. Indeed, equations (172) of appendix D shows that C-, a < 0 for a max = 90° and C 2q = 1 
for cr max = 180°. In these cases from equation (66) the pitch sensor uncertainty curve should 
concave downward for all N over from —90° to 90°. 

The results illustrated in figure 2 are summarized in columns 2 to 5 of table 1. Row 3 indicates 
the pitch angle calibration range, row 4 contains the mean normalized standard deviation over 
this range, row 5 indicates the reduced “usage range” over which measurements are to be made, 
and the final row contains the mean normalized standard deviation over the reduced “usage 
range.” Note that calibration over —45° to 45° slightly reduces the mean normalized standard 
deviation V 1 ^ 2 within the usage range over —30° to 30° compared with calibration over —30° 
to 30°. However, calibration over —90° to 90° worsens V 1 ^ 2 by 12 percent within the usage range 
from —30° to 30° compared with calibration over —30° to 30°. For calibration over —45° to 
45° or less, figures 2(a) and (b) demonstrate that the normalized curve shapes do not change 
significantly as N varies from 5 to 65. The results of figure 2 suggest that the AOA sensor should 
be calibrated over —45° to 45° degrees for use in the normal —30° to 30° range. 

The effects of unequally spaced pitch angle points within design D 0 are illustrated in figures 3 
and 4 and in columns 6 to 8 of table 1. Each calibration is conducted over a pitch range from 
—30° to 30° with 5.63° roll increments, M = 64, and N = 33. In top plot of figure 3 pitch angle 
calibration points, shown as circles, are closely spaced at 1° increments within a range from —10° 
to 10° and are more widely spaced at 4° increments for |a| > 14°. In bottom plot of figure 3, 
pitch angle calibration points are closely spaced at 1° increments for |a| > 20° and are more 
widely spaced with 4° increments for lo- 1 < 16°. Note that the normalized standard deviation 
curve of bottom plot of figure 3 is significantly flattened, although the minimum value is greater 
when compared with top plot of figure 3. Table 1 indicates that the design of bottom plot of 
figure 3 reduces V 1 ' 2 by 10 percent compared with that of figure 3 over a usage range from —30° 
to 30°; however, the latter design increases V 1 ! 2 by only 1 percent over a usage range from —10° 
to 10°. The design of bottom plot of figure 3 reduces V x !' 2 by 9 percent over a usage range of 
— 10° to 10° compared with design Dj of figure 2(a). 


13 



Figure 4 illustrates a design wherein all calibration points are located at ±30° boundaries 
except for a single center point at 0°; there is less variation of normalized standard deviation over 
the calibration interval compared with figure 3. As discussed in reference 5, designs containing 
a preponderance of boundary points reduce overall precision uncertainty at the expense of 
increased bias uncertainty due to modeling error. 

The results of figures 2 to 4 show that only small uncertainty reductions result from the use of 
nonuniformly spaced pitch calibration sets compared with design Dj. If minimum uncertainty is 
required over —10° to 10° the design of top plot of figure 3 provides a modest 9-percent average 
uncertainty reduction compared with design Dj. 

Figure 5 illustrates pitch sensor uncertainty for a modified Dj design with N = 33 and 
M = 65, with pitch angle uniformly spaced over —30° to 30°, and roll angle uniformly spaced 
over —180° to 180° with a repeated roll point at 180°. A family of normalized standard deviation 
curves is dependent on roll angle results, although deviation is small from the corresponding 
single uncertainty curve of figure 1 with design Dj. Curves are shown for 13 uniformly spaced 
roll values ranging over —180° to 180°. This modified design, convenient for experimental use, has 
insignificant disadvantage compared with design Dj. The mean normalized standard deviation 
for this case is listed in the last column of table 1. 


6. Computation of Inferred Inputs and Confidence Intervals 

6.1. Single -Axis Sensor Without Roll 

Given observed pitch sensor output v a , the corresponding inferred pitch angle a is estimated 
by inverting equation (1) so that 


a = arcsin 


The uncertainty of a is given by 


6a = 


h, ~ /■>„ 

5a 

6v a 


+ 


5„ cos (a - <p Q ) 
Then the standard deviation of a is given by 


<ts(z) 


<T eo (z) 


S a | cos (a — <f> a ) 


(76) 


(77) 


(78) 


Figure 6 illustrates the normalized standard deviation of a versus pitch angle and shows that 
inferred pitch angle uncertainty is unbounded near the extremes, a = ±90°. 

6.2. Measurements With Roll 

Given observed model attitude sensor outputs v x , v y , and v z , the corresponding inferred 
applied pitch and roll angles, a and R, are estimated by simultaneously inverting nonlinear 
equations (2), (4), and (5) as appropriate by means of Newton-Raphson iteration or other 
iterative procedure. 


6.3. Single-Axis Sensor Package With Independent Roll Measurement 

For the single-axis pitch sensor with independently measured roll angle, inferred pitch angle 
a is computed from observed sensor output v x with equation (3) as follows: 


a = arcsin 


(c jh )/5g 


cos 2 Q x + sin 2 (R + A x ) 


2/01 4 ' sin 2 fi. 


+ arctan [tan Q x sin (R + A x )] (79) 
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Thus from equation (2) and reference 5, the uncertainty of inferred pitch angle a at known 
roll angle R is given by 


8v x — S x sin Q x cos a cos ( R T A x ) 8R 

8a = ; ; 

5* [cos Q x cos a + sin Q x sin a sin (R + Axj\ 

where 8R is the uncertainty of R. The standard deviation of a is found to be 

a _ Wy X ( z )/S x + [sin Q x cos a cos (R + A x )]' 2 o%y f 2 
| cos Oj; cos a + sin Q x sin a sin (R + A%)\ 


(80) 


(81) 


where is the variance of independently measured roll angle R. If misalignment parameter Q 
is zero, the standard deviation of 5? simplifies to the following equation: 


a a ( z ) 


^Jz ) 

S x cos a 


(82) 


The inferred pitch angle uncertainty is minimum at a = 0° and unbounded near a = ±90°. 
Normalized standard deviation curves, a a /a vx , for Q,. = 1° and A x = 90°, appear in figure 7 as 
functions of a over —90° to 90° and in figure 8 as functions of R over 0° to 180°. Figure 7 contains 
two curves with measured roll angle uncertainties of 1 times and 10 times pitch sensor uncertainty, 
respectively. For these cases, inferred pitch angle uncertainty does not vary significantly with 
roll angle. Figure 8 contains three curves with measured roll angle uncertainties of 1 times, 
10 times, and 100 times pitch sensor uncertainty, respectively. Inferred pitch angle uncertainty 
varies significantly with roll only for the latter case. Note that the inferred pitch angle uncertainty 
is approximately 15 percent greater at a = 30° than at a = 0°. 


6.4. Two-Axis Sensor Package 

The two-axis model attitude sensor package containing accelerometers aligned with the x- 
and y - axes is suitable for simultaneous pitch and roll measurement within limits. As is shown, 
measurement singularities exist at ±90° pitch and near ±90° roll. Let z denote the 1x2 vector 
of inferred inputs corresponding to 1 x 2 observed output vector v, obtained by simultaneous 
solution of equations (2) and (4), where 


z = [a 
v = [v x 



(83) 


In addition, let f ( C,z) denote the 1x2 vector of functions defined by transducer equations (14) 
and (15) as follows: 

f(C,z) = [f x (c x , z) / y (c y ,z)] (84) 

where 4x2 parameter matrix C is defined as 


C = 



(85) 


The 2 x 2 J acobian matrix of equation (84) with respect to input vector z is given by 


F z 


'0f(C,z)' 


fx a 

fm' 

dz 


JxR 

fm. 


( 86 ) 
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where 


fx a = S'* [cos Fl x cos + sin Q x sin a j, sin (74, ± A a ,)] 
f XR = -5* sin Q x cos a* cos (74 + A„) 


(87) 


and 


/ y = SCcos Q v sin a* sin 74 + sin Q y (sin A y cos a ^ + cos A y sin aq, cos 74)] 
f y = ~Sy cos a^(cos Q y cos 74 — sin Q y cos A y sin R /.) 

A solution to f(C,z) = 0 exists only if Jacobian matrix F., is nonsingular at z (ref. 5). The 
singularity loci of matrix F z are obtained by setting the determinant of equation (86) to zero. 
Note that F z is singular at a = ±90°. Figures 9(a) and (b) show the singularity loci as functions 
of a and R where coning angles Q x = Q (/ equal 0.1° and 1°, respectively, for A x = 0° and 
A y = 90°; these loci nearly coincide with a = ±90° and R = ±90° for |fi| < 1°. Figures 9(c) 
and (d) illustrate the previous case repeated for A x = 90° and A y = 90°. Note the significant 
departure from R = ±90° as a approaches ±90° for |fl| > 0.1°. Parametric studies show that 
the singularity loci are dependent upon A x and nearly independent of A y for |fi| < 1°. Figure 9 
illustrates the extreme cases. 

As shown in reference 5, the uncertainty 6 z of inferred input vector z, corresponding to 
observed output vector v, is obtained from the following equation: 

<5z = 6v Ff 1 (89) 

where 4z = [8a 8R], and 8v = [8v x 5u v ] is the uncertainty of predicted output vector v. Thus 
the 2 x 2 covariance matrix of z is given by 

S z = F./'v.F . 1 (90) 

Matrix E- is the 2x2 covariance matrix of v, whose diagonal elements <r| and cr| are 

estimated by means of equation (52). Confidence and prediction intervals for z are obtained 
from equation (90). 

The normalized standard deviations of a and R, shown as <j-/ a, ;x and a j; / a vy , are presented 
for comparison in figure 10 as functions of R for selected x- and y-axis sensor output uncertainties 
as R varies from —180° to 180° at pitch angles of 0°, 20°, 40°, 60°, and 80° and at coning angles 
of 0.1° and 1°. Sensor x and y outputs are assumed to be uncorrelated; hence, is diagonal. 
As seen in the figures, inferred roll angle is singular near R = ±90°. Consequently, x-axis 
sensor misalignment correction accuracy is limited in this region, causing inferred pitch angle 
uncertainty to increase sharply near R = ±90°, although the maximum pitch error is bounded 
by coning angles Q. Roll certainties reach minima near R = 0° and 180°. 

In figure 10(a), a:-axis sensor output uncertainty equals y-axis sensor output uncertainty, that 
is, cr x = ijy] however, a y = lOcr^, in figures 10(b), (c), and (d). The a:-axis sensor azimuth A x = 90° 
in figures 10(a), (b), and (c); A x — 0° in figure 10(d). Comparison of figures 10(a), (b), and (c) 
shows that, for < 1° and |a| < 60° or 120° < !«! < 240°, the ten times less accurate y-axis 
sensor does not significantly worsen inferred pitch angle uncertainty in the ranges |7?| < 85° and 
95° < |72| < 265°. However comparison of figures 10(a), (b), and (c) shows that the inferred 
pitch angle uncertainty singularity near 90° widens as coning angle increases from 0.1° to 1° for 
(T y — 10(7, „. Figures 10(b) and (d) show that pitch angle uncertainty is least affected by roll for 
Ax = 0 °. 
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The x-y package is suitable for pitch-roll measurement in the range [|a| < 80° or 100° < \a\ < 
260°] and [|f?| < 60° or 120° < |f?| < 240°]. Note that and and or are logical operators in the 
above statement. 


6.5. Two- Axis Sensor Package 

The two-axis model attitude sensor package containing accelerometers aligned with the 
x- and z-axes is not suitable for simultaneous pitch and roll measurement at typical wind tunnel 
model test attitudes, since singularities exist near roll of 0° and ±180°, as well as at pitch of 
±90°, as shown later. Let z denote the 1x2 vector of inferred inputs corresponding to 1 x 2 
observed output vector v, obtained by simultaneous solution of equations (2) and (5), where 

z — [a ii] 

v = [v x v j 

In addition, let f (C ,z) denote the 1x2 vector of functions defined by transducer equations (14) 
and (16) as follows: 

f(C,z) = [f x ( c x ,z) /,( c.,z)] (92) 

where 4x2 parameter matrix C is defined as 


(91) 


C = 



(93) 


The 2x2 Jacobian matrix of equation (92) with respect to input vector z is given by 


’<9f(C,z )" 


’f*a 

f . ' 

dz 


. f X R 

f~n. 


(94) 


where 


f Za = 5 .{cos Qj, sin a k cos R k + sin Q z (cos A, cos a k + sin A z sin a k sin R k )] 
f ZR = S z cos a k (cos fi, sin R k — sin Q , sin A., cos R k ) 


(95) 


Figures 11(a) and (b) show the singularity loci of matrix F z as functions of a and R for 
A x = 90° and A,. = 0°, where = Q.. ranges from 0.1° to 1°; the singularity loci nearly coincide 
with the lines a = ±90°, and the lines R = 0° and R = 180° for |Q| < 1°. Figures 11(c) and (d) 
illustrate the previous case repeated for A x = 0° and A z = 0°; note the significant departure from 
R — ±90° as a approaches ±90°, for |Q| > 0.1°. Parametric studies show that the singularity 
loci are dependent upon A x and nearly independent of A z for |Q| < 1°. 

The x-z package is useful for pitch measurement from a = —180° to 180° with independently 
measured roll R except for the points {a,R} = {±90°, ± 90°}, as is now shown. Given observed 
package output v at known roll R , a is estimated by least-squares solution of overdetermined 
system (eq. (92)), where the uncertainty of the estimate is 

6S = 6vfJ(i;C) ' (96) 


and where 


3f(C,z) 

da 


[fza fz, 


It is readily shown for = 0° that 


(97) 


f Q fJ = ( 1 — cos 2 R) cos 2 a + cos 2 R 


(98) 
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for which case the estimated pitch angle uncertainty is unbounded only at the points 

{a,f?} = {i90°,i90°}. 

It is seen that the x-z package is satisfactory for pitch measurement from a = —180° to 180°, 
where roll R is measured independently, except for the points {a,R} = {±90°, ±90°}. Although 
it is capable of simultaneous pitch-roll measurement, the usable range, limited to [|«| < 80° or 
100° < |a| < 260°] and [30° < [f?[ < 150°], excludes typical wind tunnel model attitudes. 

6.6. Tin- ee-Axis Sensor Package 

The three-axis sensor package, with accelerometers aligned with the x-, y -, and z-axes, 
is suitable for simultaneous pitch-roll measurement at all attitudes, except a = ±90° where 
R cannot be determined, as shown subsequently. Let z denote the 1x2 vector of inferred 
inputs corresponding to 1 x 3 observed output vector v, estimated by least-squares solution of 
overdetermined equation system (eqs. (2), (4), and (5)), where 

z = [S fl] 1 

v = [v x v y v z } ) 

In addition, let f (C,z) denote the 1x3 vector offunctions defined by transducer equations (14) 
to ( 16) as follows: 

f(C,z) = [f x (c x ,z) / v (c v ,z)/,(c*,z)] (100) 

where 4x3 parameter matrix C is defined as 


(99) 


C = 


'c '<• 

i '- .</ 1 ■ 


(101) 


The 2x3 Jacobian matrix of equation (100) with respect to input vector z is given by 


<9f(C,z ) 


Jxa 

fib 

f*> 

dz 


JxR 

fm 

f ZR . 


(102) 


where the elements of F. are defined in equations (87), (88), and (95). 

A leasbsquares estimated solution to the 3x1 system f(C,z) = 0 exists only if ¥ z has rank 2, 
or equivalently, if 2 x 2 moment matrix F F is nonsingular. Clearly, F*Fj is singular for 
a — ±90°. General analytic computation of the remaining zeros of det (F z Fj) is unmanageable. 
However, parametric computations show that F,Fj is nonsingular for all values of R, A, and 
a ^ ±90° whenever |Q| < 10°. The singularity locus of F,F T for Q x = Q„ = Q z = 45° and 
A x = A y = A z = 90° is shown in figure 12; this case is primarily of academic interest since 
typically |Q| < 1°. 

It is shown in reference 5 that the uncertainty 5z of inferred input vector z, relative to 
observed output vector v, is obtained from the following equation as 

8z = <5vFj (FzFJ) -1 (103) 

where Sz = [6S Note that 6v = 8v y 6v z ] is the uncertainty of predicted output vector 

v. It follows that the 2x2 covariance matrix of z is given by 

S z = (F F ) 1 F,,S v Fj (FF‘) 1 (104) 


where is the 3x3 covariance matrix of v. Confidence and prediction intervals for z are 
obtained from equation (104). 
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To determine closed-form variance functions of inferred inputs a and R for the three-axis 
sensor without misalignment errors, evaluate gradient matrix F. using the parameter values of 
equation (58) as follows: 


F 

Z 


cos a sin R sin a cos R sin a 

0 —cos R cos a sin R cos a 


(105) 


Moment matrix F^Fj is then given by 


Let the y- and z-axis sensors have equal measurement variance u y and let the three 
measurement errors be uncorrelated; then measurement covariance matrix Ev is of the form 



F.F 1 


£ v = 


o 0- 

o <x 2 v o 

-0 0 <x 2 _ 


(107) 


Combine equations (104) to (107) to obtain variance functions <7?(z) and <r'|(z) of the inferred 


inputs, as follows: 


crj = <r 2 cos"’ a + a 1 sin 2 a 


,2 _ _ZJL 


(108) 


cos‘ a 


Note from equations (108) that <r~(z) = a x whenever a x = a y . If cr y > cr x then cr~(z) 
reaches a minimum of a x at a = 0°, and reaches a maximum of cr y at a = ±90°. Thus, the 
three-axis sensor eliminates inferred pitch angle uncertainty singularities at a = ±90° seen for 
the single- axis sensor with independently measured roll in equation (81) and for the two-axis 
x-y sensor package. However, inferred roll angle is unbounded at a = ±90°. Both uncertainties 
are independent of roll. 


Curves of relative standard deviations (r a (z) / a, ;x (z) and a E (z) / cr vy (z) appear in figures 12 
to 20 as a varies from —90° to 90°, as R varies from 0° to 180°, and for cr y = oq. Weighted 
least-squares estimation is assumed, where output component squared errors are weighted by 
the inverse of the associated output variances. Figures 13 and 14 illustrate inferred pitch 
and roll angle uncertainties plotted versus pitch and roll, respectively, for a y = a z = 10cr x , 
Q x = Q y = Q z = 0.1°, A x = tt/ 2 , and A y — A z — 0°. There is negligible deviation from the 
misalignment-free curves of equations (108). 


Figures 15 and 16 repeat the case of figures 13 and 14 with cr y = <j z = lOoq, except that 
Vl x = 0 ?/ = Q z = 1°; there is insignificant change from figures 13 and 14. Figures 17 and 18 
repeat the case of figures 15 and 16 with fl* = Q v = fL = 1° except that <x v = a z = a x ] 
inferred pitch uncertainty is nearly constant over pitch and roll in spite of 1° misalignment 
angles. Figures 19 and 20 repeat the case of figures 15 and 16 with cr y = cr z =. lOoq,, except that 
Q x = Q y = Q z = 5°; pitch angle uncertainty worsens by approximately 50 percent at roll of 90°. 


6.7. Summary of Pitch Measurement With Roll 

Comparison of figures 13 to 20 with figures 10 and 11 confirms that the three-axis sensor 
package is required for general purpose pitch-roll measurement. To obtain the most accurate 
pitch measurement over the full pitch and roll angle ranges, high-precision sensors are required 
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on all three axes. However, inferred pitch angle measurement accuracy can be maintained within 
the typical ranges of —60° to 60° for pitch and —180° to 180° for roll angles using y- and z-axis 
sensors whose uncertainties are up to 10 times greater than the i’-axis sensor uncertainty, and 
with sensor misalignment angles as large as 2°. Thus, accurate pitch measurement with roll can 
be obtained from —60° to 60° with a high-precision a:-axis sensor in two- and three-axis packages 
with significantly less accurate y- and z-axis sensors and in a single-axis package with significantly 
less accurate independent roll measurement. Note that roll measurement at a = ±90° is not 
possible with the three-axis sensor. The x-y axis sensor is useful primarily for pitch measurement 
from —180° to 180° with independently measured roll for R ^ ±90°. 


7. Fractional Experimental Designs 

Fractional experimental designs constructed as subsets of larger type D experimental designs 
can provide more efficient calibration while maintaining adequate prediction uncertainties. Test 
point placement for fractional designs includes the following considerations: 

1. Comprehensive test point coverage throughout the area of fi including boundaries 

2. Sufficient incremental resolution to define functional variation 

3. Limited number of experimental design points to maintain affordable calibration 


The number of points for experimental design D can be reduced while maintaining coverage 
over its full area by decimation of selected interior rows and columns. This procedure also 
maintains full incremental resolution within the nondecimated rows and columns. Figure 21(a) 
illustrates an N x M type D design, where N = 19 and M = 13. Figure 21(b) illustrates the 
same design wherein every A' fl th row is decimated by a factor of K a = 3, and every K a th column 
is decimated by a factor of K R = 4. Boundaries are not decimated. The number of points, 
denoted by C' , of the fractional design is thereby reduced from C = NM = 247 to C = 139, 
where 


a = n 


M - 1 

Kr 


+ 1 + (M — 1 ) 


1 k r 


N - 1 

Kr, 


+ 1 


(109) 


8. Replicated Calibration 

As discussed in reference 5, up to 10 replicated calibrations over an extended time period 
are necessary to obtain adequate statistical sampling over time, to estimate bias and precision 
uncertainties, and to test for nonstationarity and drift of the estimated parameters. The following 
analysis of variance techniques developed in reference 5 are applied to experimental calibration 
data presented below: 

1. Test of significance for presence of bias uncertainty 

2. Estimated bias and precision uncertainties 

3. Tests of significance for estimated offset and sensitivity drift 
Typically six replicated calibrations are obtained. 

9. Experimental Calibration Data 

Calibration residual plots are shown figures 22 to 54 for the experimental calibration data 
sets described in this section, with 95 percent calibration confidence intervals indicated as dotted 
curves and 95 percent prediction intervals indicated as dash-dotted curves. Residual sets for each 
replication are indicated by a unique symbol. Numerical statistics for selected figures are listed 
in table 2 as follows. The standard error of the regression is denoted by a E . Analyses of variance 
(ref. 5) provide estimates of standard error <7bi as due to calibration bias error and standard error 
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<r prec due to calibration measurement precision error. Symbol T bias denotes the test value for the 
calibration bias error test of significance; (-F\,ias )95 denotes the corresponding F-distributed limit 
at 95 percent confidence level. In addition, standard errors and tests of significance are indicated 
for variation between replications of estimated sensor offset and sensitivity. Variables Uj and cr<j 
denote the estimated standard errors due to drift in b and S , respectively. Symbols Tb and Ts 
denote test values for significant offset drift and sensitivity drift, respectively; (Fj.s');® denotes 
the corresponding F-distributed limit for both test values. Note that the tests are statistically 
significant if test value T exceeds limit F. 

Inferred residual plots are also provided for each data set, obtained by back-computation of 
inferred calibration inputs using the observed calibration output data and estimated calibration 
parameters. The corresponding inferred calibration confidence intervals and prediction intervals 
are shown as dotted curves and dash-dotted curves, respectively. 

9.1. Single-Axis Calibration Without Roll 

Figures 22, 24, and 26 present calibration residual plots with 95 percent confidence and 
prediction intervals for six replicated calibrations without roll of a high-precision single-axis 
AOA sensor, without temperature correction. Inferred calibration inputs are back- calculated by 
using equation (76). The corresponding inferred residual plots appear in figures 23, 25, and 27. 

The calibration of figures 22 and 23 employs design Dj from a = —36° to 36° with 
2° increments. The standard error of regression of figure 22, listed in table 2, is 0.000160°; no 
significant calibration bias error or sensor sensitivity drift over the six replications is detected. 
Slightly significant sensor offset drift is detected. The rms value of the residuals of the inferred 
angles, denoted by <r inv , equals 0.000174°. 

The calibration of figures 24 and 25 employs design Dj from a = —180° to 180° with 
5° increments. The calibration residuals disclose a systematic sinusoidal error pattern with 
two periods from a = —180° to 180°. Note in table 2 that the larger standard error of regression 
for figure 24 is 0.000317°, compared with figure 22, and significant calibration bias uncertainty is 
detected. Significant sensor offset and sensitivity drift are not detected. At ±90° where inferred 
confidence and prediction intervals become unbounded, most residuals of the inferred angles fall 
outside the boundaries of figure 25. The observed sinusoidal systematic error in figure 24 is due 
to static deflection of isolation pads within the sensor package. 

Figures 26 and 27 illustrate residuals for six replicated four-point tumble calibrations with 
a standard error of regression for figure 36 of 0.000284° listed in table 2. The large calibration 
confidence intervals are caused by the reduced number of degrees of freedom. Note also that 
significant calibration bias uncertainty is detected although without significant sensor parameter 
drift. Most residuals of the inferred angles fall outside the chart boundaries at a = ±90° in 
figure 27. 

9.2. Single -Axis Calibration With Roll 

Two single-axis AOA sensors were simultaneously calibrated with roll over multiple repli- 
cations. Sensor 1 is a high-precision unit; sensor 2 is a less expensive unit of lower accuracy. 
Experimental design Dj with an extra roll point at 180°, as in the design of figure 5, was employed 
with pitch angle limits of ±30° and ±180°. 

9.2.1. Full calibration from. —30° to 30°. Pitch and roll angle step sizes are 5° and 15°, 
respectively, and the resultant design contains 325 calibration points per replication over six 
calibrations. Temperature variation did not exceed 1°C during calibration. 

Figure 28(a) illustrates calibration residuals of sensor 1 computed without temperature 
correction; residuals are plotted versus pitch angle. As seen in table 2, the standard error of 
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the regression is 0.000776° with only minimally significant indicated calibration bias uncertainty. 
However, very significant sensor sensitivity drift, with T 9 = 919, and less significant sensor offset 
drift are detected, which is also apparent from slope variations seen in the residual pattern. 

Figures 28(b) and 29 illustrate calibration residuals and residuals of the inferred angles, 
respectively, for the dataof figure 28(a) recomputed with temperature corrections for sensor offset 
and sensitivity. Standard error reduces to 0 .000387° compared with that in figure 28(a), as shown 
in table 2; significant calibration bias uncertainty is detected. After temperature correction, 
sensor offset drift and sensitivity drift are greatly reduced, with T 9 = 5.14. Figure 30 illustrates 
individual residual curves for the first replication only plotted versus pitch and parameterized 
by calibration roll angles from 180° to 0 ° by using calibration parameters estimated over six 
replications. The systematic error pattern produces minimum error dispersion at —5° pitch and 
greatest dispersion at ±30° pitch. Figures 31 and 32 illustrate calibration residuals and residuals 
of the inferred angles, respectively, plotted versus roll angle. Minimum dispersion is apparent 
near ±90° roll, with maximum dispersion near 0° and ±180° roll. Statistics for figure 31 are 
identical to those for figure 28(b). Figure 33 illustrates individual inferred residual curves for 
the first replication only plotted versus roll angle and parameterized by calibration pitch angles 
over —30° to 30°, using parameters estimated over six replications. 

Figure 34 illustrates calibration residuals for less accurate sensor 2 plotted versus pitch with 
temperature correction. The standard error of the regression is 0.00166° as listed in table 2; 
calibration bias uncertainty is insignificant. Strongly significant sensor offset and sensitivity 
drifts are indicated, which are apparent in the residual patterns. 

9.2.2. Fractional calibration from —30° to 30°. The design cardinality of the 325-point 
calibration D design in section 9.2.1 is reduced to 53 points as follows: overall pitch and roll 
angle resolutions are reduced from 3° to 15° and from 15° to 30°, respectively. Alternate rows 
and columns are then decimated by factors of 2. Figure 35 illustrates the fractional calibration 
residuals for sensor 1. Note the enlarged calibration confidence intervals, caused by reduced 
degrees of freedom, and the larger prediction intervals compared with the full calibration data 
of figure 28(b). As seen in table 2, the standard error is increased from 0.000387° in figure 28(b) 
to 0.000427°; the test for calibration bias error is significant. 

Figure 36 illustrates the data residuals computed from the full data set by using parameter 
vector c and confidence intervals obtained from the fractional calibration. The standard error 
of the residuals equals 0.000389° compared with the standard error of 0.000387° obtained for 
the full data set of figure 28(b). For sensor 1, calibration by this particular fractional design 
provides a fit nearly equivalent to that provided by the complete design. 

9.2.3. Calibration from —180° to 180°. Pitch and roll angle step sizes are f5° and 
30°, respectively, with 325 calibration points per replication. Temperature variation during 
calibration did not exceed 1°C. 

Figures 37 and 38 illustrate sensor 1 calibration residuals and residuals of the inferred angles, 
respectively, computed with temperature correction over four replications; residuals are plotted 
against pitch angle. The standard error of the regression is 0.000489° with significant indicated 
calibration bias uncertainty, as seen in table 2. Slightly significant sensor offset drift is detected 
without significant sensor sensitivity drift. Figure 39 illustrates individual residual curves for 
the first replication only using calibration parameters estimated over four replications; curves 
are plotted versus pitch angle and parameterized by calibration roll angles from 0° to 180°. The 
systematic residual pattern is dependent on both pitch and roll; error variation with pitch angle 
is sinusoidal with two periods over a = —180° to 180°. 

Figure 40 illustrates calibration residuals plotted versus roll angle. Figure 41 illustrates 
individual residual curves for the first replication using calibration parameters estimated over 
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four replications; curves are plotted versus roll angle and parameterized by calibration pitch 
angles from 0° to 180°. The systematic error pattern is dependent on both pitch and roll; error 
variation with roll angle is sinusoidal with one period over R = —180° to 180°. 

Figures 42 and 43 illustrate sensor 2 calibration residuals and residuals of the inferred angles, 
respectively, with temperature correction over six replications; residuals are plotted versus pitch 
angle. The standard error of the regression is 0.00134°. Other statistics appear in table 2. 
Calibration bias uncertainty is insignificant. Strongly significant sensor offset and sensitivity 
drift are detected between replications. 

9.3. Tliree-Axis Calibration With Roll 

A three-axis model attitude sensor package containing identical high-precision sensors was 
calibrated with roll for six replications. Experimental design Di with an extra roll point at 180°, 
as in the design of figure 5, was employed with pitch angle limits of ±90° and ±180°. Sensor data 
are temperature corrected; confidence and prediction intervals appear in each figure. Residuals 
of the inferred angles are obtained by subtracting true angle values from the back-computed 
angle values. 

9.3.1. Calibration from —90° to 90°. Pitch and roll angle step sizes are 10° and 30°, 
respectively, with 247 calibration points per replication over six calibrations. Total calibration 
time was approximately 13 hr with temperature variation no greater than ±1°C. Figures 44(a), 
(b), and (c) illustrate calibration residuals plotted versus pitch angle for the x-, y-, and z-axis 
sensors, respectively, over six replications. The regression standard errors of the three sensors 
are 0.000434°, 0.000444°, and 0.000355°, respectively. As seen in table 2 significant calibration 
bias uncertainty and significant offset drift are detected for each of the three sensors. However, 
significant sensitivity drift is detected only for the x- and a- ax is sensors. 

Figures 45 and 46 illustrate residuals of the inferred pitch and roll angles, respectively, for 
the first replication only; curves are plotted versus pitch angle. Prediction intervals for inferred 
roll angle uncertainty, shown as functions of pitch angle, are significantly greater than those for 
inferred pitch angle uncertainty. 

9.3.2. Calibration from —180° to 180°. Step sizes for pitch and roll angles are 10° 
and 30°, respectively, with six replications. Total calibration time was approximately 28 hr. 
Figures 47(a), (b), and (c) illustrate calibration residuals for the x-, y-, and z-axis sensors, 
respectively, over the six replications; curves are plotted versus pitch angle. Statistics are given 
in table 2. The regression standard errors of the three sensors are 0.000409°, 0.000523°, and 
0.000479°, respectively. Significant calibration bias uncertainty is detected for each of the three 
sensors. Two periods of a sinusoidal error pattern over a = —180° to 180° are apparent in 
figure 44(c) for the z-axis sensor. However, significant sensitivity drift and significant offset drift 
are detected only for the x- and z-axis sensors. Figures 48 and 49 illustrate inferred pitch and roll 
angle residuals, respectively, for the first replication only; curves are plotted versus pitch angle. 
A sinusoidal error pattern is also apparent in the inferred pitch angle residuals of figure 48, with 
unusually large scatter at a = —90°. The observed sinusoidal systematic error is due to static 
deflection of isolation pads within the sensor package as observed also in figure 24. 

9.3.3. Six-point tumble calibration. Six replicated six-point tumble calibrations using 
design T are obtained from the previous data set in section 9.3.2. Figure 50 illustrates x-, 
y-, and z-axis sensor calibration residual curves over the six replications; individual curves are 
plotted versus pitch angle. Statistics appear in table 2. The regression standard errors are 
0.000433°, 0.000197°, and 0.000279°, respectively. Significant calibration bias uncertainty is 
detected for the x- and y-axis sensors. However, neither significant offset drift nor sensitivity 
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drift is detected for the x- and y-axis sensors. Figure 51 illustrates inferred pitch angle residuals 
and roll angle residuals, respectively, for all six replications; curves are plotted versus pitch angle. 

Figure 52 illustrates sensor output residuals for the entire calibration data set computed 
by using parameters estimated from the six-point tumble calibration data. The indicated 
confidence and prediction intervals are obtained from the tumble calibration regression analysis. 
Standard residual errors are 0.00104°, 0.00095°, and 0.00084°, respectively. The corresponding 
regression standard errors appear in the previous paragraph. Comparison with figure 44 shows 
that the replicated six-point tumble test significantly underestimates prediction intervals. At 
the same time it suffers greater calibration uncertainty compared with the full calibration, as 
evidenced by the larger calibration confidence intervals. Compared with figure 47(a), figure 52(a) 
illustrates increased standard residual error (0.00104° compared with 0.000433°), as indicated 
by the systematic error pattern, caused primarily by the limited spatial resolution of the T 
experimental design compared with the multipoint Do design. 

9.3.4 ■ Fractional calibration from —180° to 180°. The design cardinality of the 481-point 
calibration of section 9.3.4 is reduced to 73 points as follows: overall pitch and roll angle 
resolutions are reduced from 10° to 30° and from 30° to 60°, respectively. Alternate rows and 
columns are then decimated each by a factor of 2. Statistics are given in table 2. Comparison 
of the x-axis sensor fractional calibration residuals, shown in figure 53, with the full calibration 
residuals of figure 47(a) shows nearly the same prediction intervals, although the calibration 
confidence intervals are enlarged due to fewer degrees of freedom. The standard errors and 
tests for calibration bias error are nearly unchanged. However, the fractional calibration fails 
to detect significant offset drift and indicates considerably reduced sensitivity drift significance. 
Figure 54 illustrates the data residuals computed from the full data set by using the parameter 
vector c and confidence intervals estimated by fractional calibration. The standard error of 
the residuals shown in figure 54 equals 0.000410° compared with a residual standard error of 
0.000409° obtained in figure 47(a). Except for offset drift detection, the 73-point fractional 
calibration performs equivalently to the full 481-point calibration. 

10. Concluding Remarks 

Statistical tools, developed in NASA/TP-1999-209545 for nonlinear least-squares estimation 
of multivariate sensor calibration parameters and the associated calibration uncertainty anal- 
ysis, have been applied to single- and multiple-axis inertial model attitude sensors with and 
without roll. These techniques provide confidence and prediction intervals of calibrated sensor 
uncertainty as functions of applied input angle values. They also provide a comparative per- 
formance study of various experimental designs for inertial sensor calibration. The importance 
of replicated calibrations over extended time periods has been emphasized; replication provides 
estimates of calibration precision and bias uncertainties, statistical tests for calibration or mod- 
eling bias uncertainty, and statistical tests for sensor offset and sensitivity drift during replicated 
calibrations. 

The techniques developed herein properly account for correlation among estimated calibration 
parameters and among multisensor signal conditioning channels, allow inclusion of calibration 
standard uncertainties, and account for uncertainty of independently measured roll angle. Previ- 
ous empirical techniques for treating correlations among estimated parameters may overestimate, 
or in certain cases significantly underestimate, uncertainty magnitudes. 

The sensor output variance function, and hence calibration confidence intervals and prediction 
intervals, have been shown to be identical for x-, y-, and z-axis sensors. Moreover, the output 
variance function is independent of the inertial sensor parameters c = [feSTIA] 1 . Hence, the 
design figure of merit is independent of the sensor under calibration. In addition, the sensor 
output variance function is independent of roll angle R for experimental design Do, wherein roll 
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angle test points are uniformly spaced over the roll angle range without repeated principal angle 
values. 

Parametric studies show that the pitch sensor figure of merit, computed within a limited usage 
range, can be reduced by limiting pitch angle test points to a range approximately 1.5 times 
the usage range. For example, calibration over a pitch range from —45° to 45° is appropriate 
for a pitch usage range of —30° to 30°. Additional modest variance reduction within a limited 
test range is possible by concentrating pitch angle test points near the center of the range of 
interest. However, as discussed in NASA/TP-1999-209545, uniformly spaced designs minimize 
the mean normalized error variance due to systematic bias errors. For this reason, design Dj 
with uniformly spaced pitch and roll angle test points is preferable. Experimental results show 
that calibration over a pitch range from —180° to 180° detects systematic bias errors not seen 
in pitch calibrations from —45° to 45°. 

Experimental results show that fractional multipoint D designs can provide adequate statisti- 
cal uncertainty and uncertainty characterization with increased calibration efficiency. However, 
experimental results show that tumble test T calibration designs, limited to cardinal angles, pro- 
vide insufficient spatial resolution to adequately characterize systematic modeling uncertainty. 
As a result, prediction intervals tended to be significantly underestimated in spite of increased 
calibration uncertainty due to fewer degrees offreedomevidenced by larger calibration confidence 
intervals. 

Simple closed-form rational trigonometric polynomial expressions are obtained for computa- 
tion of confidence and prediction intervals for design Dj . In any case, numerical poinhby-point 
calculation of confidence and prediction intervals for any design is readily programmed for on-line 
computation or posttest data reduction. 

Inferred input pitch and roll angle uncertainties are dependent upon independent variables, 
pitch angle a and roll angle i?, for any experimental design, even if the variance function is 
independent of R. 

Single- and two-axis model attitude sensors do not provide accurate pitch angle or roll 
angle measurements near pitch of ±90°. Neither does the two-axis sensor provide accurate roll 
measurement near roll of ±90° at any pitch angle. Within the range of typical sensor parameters 
the three-axis sensor eliminates measurement singularities except for roll angle measurement 
near pitch of ±90°. By using identical x-, y -, and z-axis sensors, full pitch angle precision is 
maintained over a pitch range from — 180° to 180°. Adequate pitch angle measurement precision 
with roll can be maintained within a pitch angle range from —60° to 60° by use of a precision 
X-axis sensor with significantly less accurate y- and z-axis sensors, such as a y = u z — 10er„, and 
O < 2°. 

Recommendations for model attitude sensor calibration and usage are as follows: 

1. The pitch angle calibration range should be approximately 150 percent of the usage range. 

2. The roll angle calibration range should be from —180° to 180°. 

3. Test points should be uniformly spaced in both pitch and roll. 

4. Pitch angle should vary from minimum angle to maximum angle and back to minimum 
angle. 

5. Fractional D calibration experimental designs may be employed for calibration efficiency, 
provided that statistical adequacy is established experimentally. 

6. Calibrations should be replicated at least 6 times, and preferably 10 times, for estimation 
of bias and precision uncertainty and for detection of parameter nonstationarity. 
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7. Four-point and six-point tumble calibration experimental designs are not recommended 
for laboratory calibration. 

8. The single-axis package may be used for pitch angle measurement with adequate uncer- 
tainty whenever the uncertainty of the independently measured roll angle does not exceed 
10 times the desired pitch angle uncertainty. 

9. The three-axis sensor package is suitable for general pitch-roll measurement with adequate 
accuracy except for roll measurement near pitch of ±90°. The y- and z-a.xis sensor 
uncertainties should not exceed 10 times the a;- ax is sensor uncertainty. 

10. The x-y axis sensor package is suitable only for measurements away from pitch of ±90° 
and roll of ±90°. The y-axis sensor uncertainties should not exceed 10 times the x-axis 
sensor uncertainty. 

The recommended calibration experimental designs may be readily implemented by means of 
modern automated calibration apparatus. 
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Appendix A 


Derivation of y-, and z-Axis Sensor Outputs for Measurement With Roll 

The inertial attitude sensor output is obtained in reference 2 by computation of the projection 
of the gravitational force vector onto the sensor sensitive axis. The effects of package rotations 
in pitch, roll, and yaw, as well as package misalignments Q and A, are computed by means of 
coordinate transformations. 


Consider a three-dimensional right-hand coordinate system with axes x, y , and z, where 
negative z represents the direction of gravity in gravitational coordinates, shown in figure Al. 
Let x denote the direction of the model axis in model coordinates at zero pitch, roll, and yaw. 
Then g = [0 0 — 1] T denotes the normalized gravitational force vector in gravitational coordina tes, 
and let g q = [g qv g qij g q J T denote g transformed into sensor coordinates. 

Transformation from gravity coordinates to model axis coordinates, and thence to sensor 
coordinates, consists of an ordered sequence of rotations, defined by the following coordinate 
transformations: 


1. Pitch a — left-hand rotation about y-axis: 

T„(a) = 


cos a 0 —sin a 

0 1 0 

L sin a 0 cos a J 


2. Roll R — left-hand rotation about a: -ax is: 


T r (R) = 


1 0 0 

0 cos R sin R 

L 0 —sin R cos R J 


3. Yaw Y — left-hand rotation about z-axis: 


T y(Y) 


cos 1’ sin Y 0 
— sin Y cos Y 0 

0 0 1J 


(110) 


( 111 ) 


(112) 


Model Attitude Transformation 

Let the model be oriented at pitch angle a and roll angle R. Transformation from gravity 
coordinates to model coordinates is represented by pitch rotation T a (a) followed by roll rotation 
Tif(if). Gravity vector g transformed to model coordinates becomes 

gM = Tfl(i?) T a (a)g (113) 


Transformation to x-A xis Sensor Coordinates 

The sensitive axis of the x-axis sensor is nominally aligned with the model a:-axis. Sensor 
misalignment is represented as transformation from model coordinates to sensor coordinates as 
positive roll rotation T fl (A) followed by positive yaw rotation Ty(fl). In ar-sensor coordinates 
the gravity vector is given by 

gc„ = T Y (n x )T B (A x )g M (114) 

The ^-component of g^,, corrected for sensor sensitivity S x and offset b x , yields equation (2). 
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Transformation to y - Axis Sensor Coordinates 

Transformation to the sensitive axis of the y-axis sensor, nominally aligned with the model 
y-axis, is represented by the y-component of vector g m- Sensor misalignment is represented by a 
model-to-sensor coordinate transformation as positive pitch rotation T a (A) followed by positive 
roll rotation T^(f2). In y-axis sensor coordinates, the gravity vector is given by 

&to = T fl (n y )T a (A y )g Jtf (115) 

The y-component of g qj( , corrected for sensor sensitivity S y and offset b y , yields equation (4). 

Transformation to 2 - Axis Sensor Coordinates 

Transformation to the sensitive axis of the z-axis sensor, nominally aligned with the model 
z-axis, is represented by the z-component of vector g, v/ . Sensor misalignment is represented by a 
model-to-sensor coordinate transformation as positive yaw rotation Ty(A) followed by positive 
pitch rotation T a (Q). In z-axis sensor coordinates the gravity vector is given by 

g„ = 1 T v ( A. )g w (116) 

The z-component of g q _, corrected for sensor sensitivity S~ and offset b z , yields equation (5). 



Figure Al. Cartesian coordinate system. 
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Appendix B 
Evaluation of Matrix 

Matrix H £ of equation (45) is evaluated for the pitch sensor. The fcth 4x4 matrix contained 
in 4 x 4 x A' array F cc defined in equation (46) is equal to the Jacobian matrix of equation (2) 
evaluated at the fcth element of experimental design 6. The elements of are obtained for the 


pitch sensor by differentiating equations (21) to (24) as follows: 

fbb = fbS — fbQ = fbA = fss = 0 (H7) 

fsQ, = — sin Q sin a k — cos Q cos aq sin (R k + A) (118) 

fsA k — —sin £2 cos aq cos ( R k + A) (119) 

/jjjj = — 5 [(cos Q sin aq — sin Q cos aq sin (R k + A)] (120) 

foA = —S cos cos a k cos(R k + A) (121) 

f AAi = S sin Q cos a k sin (Aq + A) (122) 


Similar expressions result for the roll sensor. Matrix F cq is therefore of the form 

0 0 0 0 
0 0 fsQ k fsA k 

0 fsQ k J'VA:, f(U k 

0 fsA k fciA k f-A, 

If measurement covariance matrix U y equals <r 2 I, then matrix is given by 

0 0 0 0 
K K 

0 0 y: e k f S o j y~]et/.<?Ajr. 

k=\ k = i 

K K K 

0 ^ e hfsQ k e kfsiOt ^ e kfoA k 

k=l k=l k=\ 

K K K 

0 ^ ef/s'.q. X/'w./t h 

&=i &=i &=i 

where eq. is the fcth element of residual vector e. 

Simulation studies show that among the experimental designs evaluated above the prediction 
uncertainty is unaffected by matrix for values of measurement error standard deviation 
a < 0.01 1 v| , where v is the observed output vector. Moreover, only insignificant random 
effects are evident for a < 0 . 1 1 v | ; this confirms that H £ may be neglected for typical levels 
of measurement error. 





(123) 
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Appendix C 

Properties of Sensor Variance Functions 

The proofs of Theorems I to IV are given in this appendix. 

Theorem I: Sensor output variance function cr'lfz) is independent of calibration parameters b 
and S. 

Proof : It is shown in appendix D, in the general evaluation of elements of matrix R a . for the 
x-axis sensor with roll, that matrix It., is of the form 


R: 


1 

-J 

r bS 

r b Q 

r bA " 


nb 

ns 

SpbQ 

Swp hA 

ns 

r SS 

rsa 

rsA 


ns 

r ss 

SpsQ 

SwpsA 

rba 

rsa 

ran 

rQA 


SpbQ 

S PSQ 

S 2 pan 

S 2 WpQA 

1 

< 3 ^ 

r SA 

r QA 

r aa . 


_Swp bA 

Swp SA 

S 2 wp QA 

S 2 w' 2 p AA _ 


(125) 


The terms denoted by r and p are obtained by means of equations (130) to (135) and are 
explicitly evaluated in appendix D by equations (176) to (205). It follows from equations (28) 
to (40) for sensors y and z that matrices R„ and R may be expressed in the same form. For 
the x-, y -, and z-axis sensors, vector f c is of the form 


f c T = [fb fs fa f a ] = [fb fsS(f> n Sw4>a 


(126) 


where <fn and <f A are independent of 6 and S ; presubscripts x, y , and z are elided for convenience. 
It is shown by Lemma 1, appendix D, that if matrix R"' 1 exists, then 


t) 


where 

and 


* 9r = fc' K ' fc = P 1 </>c 

<t>c = [fb fs 4>a 4 >aY 


(127) 

(128) 

(129) 


Because the elements of (f> c , and those of P as shown, are independent of parameters b and S , 
variance function is independent of b and S . In particular, matrix P is obtained as 

P = $J$ C (130) 


r bb 

ns 

PbQ 

PbA 

ns 

rss 

PSQ 

PSA 

PbQ 

PSQ 

Pqq 

Pqa 

_ PbA 

PSA 

Pqa 

Paa _ 


where K x 4 matrix is defined as 


where 


r ^.t 


= 


K 

*1 


>.9 \ 4>n \ 4 >a 


L 4 > ’ 


c Ii ' 


4>b = [/&! • • • fb K V 


( 131 ) 


(132) 
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ll 

•■/%F 

(133) 

-e- 

£> 

II 

» 


(134) 

Z®: 

II 

-$■ 


(135) 


Therefore matrix P is independent of parameters b and S because P is computed by using 
equations (130) to (135). The proof for the single-axis sensor without roll is analogous. 

QED 


Theorem II: Sensor output variance function cr y (z) is independent of calibration parameter Q 
Proof: For x-, y -, and z-axis sensors define vectors 

gf = [ 1 sin a cos a (sin R cos A x -(- cos R sin Af) —cos o- (cos R cos A x — sin R sin A,,) ] 

g T c = [ 1 —cos a sin R —sin A y sin a + cos A y cos a cos R cos A y sin o- + sin A y cos a cos R 

g 2 T c = [ 1 —cos a cos R —cos A z cos a + sin A, cos a sin R —sin A z sin a — cos A. cos a sin R 

(136) 

and matrix 

'10 0 0 

0 cos 0 —sin Q 0 

r w = 

0 —sin Q —cos Q 0 

.0 0 0 1 

Matrix T w is orthogonal; that is, 

rl-r w = r w rl =i 


(138) 


(137) 


For sensors x, y, and z, gradient vector <f> c , defined in equation (128), equals the product of 
vector g c and matrix r as follows: 

4>1 = gjr-r (139) 

Similarly, K x 4 matrix $ c defined in equation (131) may be written as 

*c =G f r,, (140) 


where K x 4 matrix G c is defined as 


G, = 


rT 

»ci 


g, 


C K 


g b ! g.? ! go ! g.4 


(141) 


where for sensor x 


g* = [i--ir 


g, , = [sin oq • • • sin a K J 




cos cq sin ( R l + A x ) 

L cos a K sin ( li K + A x ) J 


(142) 

(143) 

(144) 
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COS OLy COS (R y + A x ) 


Kx, = 


cos a h - cos 


{Rr + As) - 


Vectors g % , g vs , g m , g V4 and g^, g z? , g ZQ , g, A are defined similarly. After noting that 

p 1 = r w- (gJ g c ) _1 r w (146) 

it follows that the gradient vector cf> c obtained in equation (128) may be combined with 
equations (139), (140), and (146) to yield 

4 * 4>? P 1 4>c = gj (Gj G e ) - 1 g c (147) 

a E 

Therefore, cr[;(z) is independent of Q for sensors x, y , and z . The proof for the single-axis sensor 
without roll is analogous. 

QED 

Theorem III: Sensor output variance function o^( z) is independent of calibration parameter A. 
Proof: Define vectors 


h„ c = [ 1 sin a cos a sin R —cos a cos R] 1 

hj, c = [ 1 —cos a sin R cos a cos R sin a j 1 

h, c = [ l —cos a cos R —sin a —cos a cos R] 1 


and matrix 


‘10 0 O' 

0 10 0 

0 0 cos A —sin A 

0 0 sin A cos A 


Note that matrix is unitary and that 

gj = lb' P-4 

for sensors x , y , and z. Define K x 4 matrix H c similarly to G c as 


After noting from equations (141) and (149) that 


it follows that 


[Gj G c ] 1 = rj [Hj H c ] 1 Ta 
4 ~ P 1 <t>c = gl [Gj G e ] -1 g c = hi [Hi H c ] 1 h e 
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for sensors x, y, and z. Therefore, af(z) is independent of azimuth A since li c and H r are 
independent of A. 

QED 


Theorem IV: Let roll angle calibration set B fl , defined in equations (10), contain K = NM points 
uniformly spaced over the interval [— 7r,7T — A 12], where M and N are integers, 
A R = 27 r/M, and the principal value of each angle contained in B/j occurs with 
the same frequency, then the pitch sensor output variance cr'^(z) is independent of 
roll angle R. 

Proof: Since variance function cP(z) is independent of calibration parameters b, S, £2, and R, 
evaluation of equations (21) to (26) using the parameter values of equations (58) yields the 
following equations: 

U = 1 

fx s = sin a 
<j> Xa — —cos a sin R 
<f> XA = —cos a cos R , 

Evaluation of equations (173) and (174) in appendix D yields 



C'mr = 0 ) 
C-2B =0 J 


It follows from equations (176) to (203) that 


where 


It then follows that 


’A w , A-,., 0 0 

r *6.9 r *.9.9 0 0 

0 0 PxQQ 0 

- 0 0 0 P*AA- 


r,„ = MN 

Hb 

r xis = M S a 


r xss = —M(N — C 2a ) 


OCfifi 


Q. 

II 

— i 1 ^ 

II 

r 

— T 

*ss 



D 

D 

~ r ^.9 

I3L 

D 

D 


0 

1 

P*AA 

0 


(155) 


(156) 


(157) 


(158) 
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where D = r Xbl r xt}S — rj Evaluate equation (59), with the help of equations (128), (154), 
and (158) to obtain 



4>c{ z)P l $(z) 

r, — 2r„.„ sin o- + r,„ sin 2 a 
D + cos 2 a/(p XAA ) 


(159) 


which is seen to be independent of roll R. 
QED 
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Appendix D 

Evaluation of the Moment Matrix 

Lemma 1: Proof of equation (127). 

Define matrix £ as 

'10 0 O' 

0 1 0 0 
0 0 S 0 
.0 0 0 Sw_ 

It follows from equations (126) and (128) that 

fc = ^4>c 

and from equations (125) and (129) that 

R = EPE 


If R 1 and £ 1 exist, then 

R 1 = £ 'P '£' 1 

Hence, 

q R = fj R J f c = P 1 £- J £ (f> c = P 1 <f> c 

QED 


(160) 

(161) 

(162) 

(163) 

(164) 


The following definitions and relations are used in the subsequent development: 
1. Pitch angle set B Q contains N points in the closed interval [cv 1Tlill cv max ] 


a. 


c. With these definitions, 


N 

Sa — ^2 sin a n 

n—1 

N 

C a — ^ ^ cos oc n 

n—1 

N \ 

s 2a = ^2 sm 2a " 

n—1 

► 

N 

C- 2a = ^ cos 2a„ 

n—1 J 


_£^ i 

^2 sin 2 a n = - C 2a ) 

n—1 

> 

N l 

^cos 2 a n = -(AT+ C- 2a ) 

n=l J 


(165) 


(166) 


(167) 
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2. Roll angle set 6^ contains M points in the closed interval [-Rmin-^max] 


a. 


M 


S R = sin R ir 

m= 1 
M 

Cr = cos R n 


b. 


m=l 


M 


S-2R = Sl " 


m=l 

M 


Cm = X^ cos 2R„ 


m= 1 


c. Using these definitions gives 

M 


J^sin 2 R m = M - C 2R ) 

m=l 

M . 1 
X^cos 2 R m = — (M + Cm) 


m = 1 


( 168 ) 


(169) 


(170) 


Special Experimental Designs 

Minimal design D 0 : Roll angle set 8^ contains M points uniformly distributed over the closed 
interval [— 7r+Ai?,7r] where the principal value of each angle contained in 6^ occurs only once 
and A R = '2i r/M . 

Minimal design Dj C D 0 : Pitch angle set B a contains N points uniformly distributed over 
the closed interval [— a m ax Q: max] where the principal value of each angle contained in fi a occurs 
only once unless cr max = 7 r; Act = 2ct max / (N — 1). Roll angle set 8^. is the same as in design D 0 . 

For D n and Dj designs containing M D copies of a minimal design, the expressions obtained 
below are multiplied by M D . 

For design Di C Dq. 


1. 


N 

Sa — sin ct„ = 
71=1 


N 

c^x: c ° s a » = 

n—l 


0 


sin[Afa imx /(AT - 1)] 
sin[a max /(AT - 1)] 


(171) 


2 . 


N 

S 2a = X^ sin 2 “n = 0 

n—l 


N 

C 2a = X^ COS 2a n 

n—l 


sin [2 IV Q max /(AT — 1)] 
sin[2a max /(Ar - 1)] 




( 172 ) 
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The following expressions are evaluated for containing M points uniformly distributed 
over the closed interval [— R^^+AR ,-R max ], where A R = 2f? max / M . Evaluation at f? max = ir for 
design D 0 yields zero in each case. 


1. 


M 


Sr = sin R m = sin = 0 


m = 1 
M 


Rr 


Cr = ^2 COS R m — cot 22 sin Rmax = 0 


m=l 


(173) 


2 . 


M 


S- 2 R — ^ ^ sm 2R m — sm 27?. max — 0 


m=l 

M—l 


2 R., 


C-2R = ^2 cos ‘ 2R »> = cot 17 sin = 0 


m = 1 


M 


(174) 


Evaluation of x-Axis Sensor Moment Matrix 

Moment matrix R, ( ,, defined in equation (44) and required for computation of variance 
function cr/(z), is now evaluated in general using the approximation in equation (49) showing 
that R r may be expressed in the form of equation ( 125), as needed for proof of Theorem I. Since 
<j';(z) is independent of parameters 6, S’, 0, and A, matrix R, is then simplified by using the 
values given in equations (58) for later evaluation of cr/(z). Further simplifications are obtained 
for designs D, D n , and/or Dj. 

With the values of equations (58), the elements of gradient vector <f> Xc become 


fn = l 


fx s = sin a 

4> xq = —cos a sin R 

4> Xi = — cos a cos R t 

Element -by- element evaluation proceeds as follows: 

1- m 

General evaluation using equations (58): 


*66 


= ffi 





(175) 


(176) 


Specific evaluation for design D: 


r x ,, = M N 

bb 


( 177 ) 


2 . ns = rsb 
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General evaluation using equations (58): 


~ ~ f x Sk 


*=i 


K 

sin a?k 
k= 1 


(178) 


Specific evaluation for design D: 


M N 

r * 6 .9 — EE sin a n — MS a 

m=l n=l 


(179) 


Specific evaluation for design Dj: r xiC . = 0. 

3- r m = 


General evaluation: 


K 

~ — ^si. — S x Pz Kt 

k=\ 


Using equations (58) gives 


(180) 


A' 



*=1 


K 

— cos ttj. sin f?.;,, 
t=i 


(181) 


Specific evaluation for design D: 


M N 

p XK) — - sin Rm cos a " = — SrCc 

m=l n— 1 


(182) 


Specific evaluation for design D 0 : 


P*ki ~ ® 


4- r iA = r M 
General evaluation: 

K 

^X‘ 0 4 — 1(‘/.1( i — ( 4* Xj^j. — S x WxPx bA (183) 

*=i 

Using equations (58) gives 


A* K 

P*bA = = - ^cos a k cos (184) 

k=l k=l 


Specific evaluation for design D: 


N M 

Px iA = - y^cos a n y^cos II, „ = -CcCr (185) 

n— 1 m — 1 

Specific evaluation for design D n : 
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(186) 

( 187 ) 

( 188 ) 

( 189 ) 

( 190 ) 

( 191 ) 

( 192 ) 

( 193 ) 

( 194 ) 
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Specific evaluation for design D 0 : 


P x sa ~ 0 


8- 

General evaluation: 

r - f f 

x aa zq *n 

Using equations (58) gives 


5? 


E 

k=l 


> 2 


= sln 




PzQSl ~ 


y 

^ x Qk 


k=l 


K 

y^^ cos 2 a k sin 2 R k 
k= 1 


Specific evaluation for design D: 

AT A' 

Pzm = X] sin2 ^ cos2 “» 

m=I n=l 


-( m - <■',„)( .v + C 2a ) 


Specific evaluation for design D 0 : 

/V,, = \m{n + c 2a ) 

9- r QA = r 4fi 
General evaluation: 


K 

T 'xQA ~ ^ ftxQ 1-4 * xai- ~ S x W x p xQA 

k=L 

Using equations (58) gives 

k i K 

P x qa — ^ ^ ik 77 ^ ~ cos* ot k sin 2R k 

k=l k = 1 

Specific evaluation for design D: 

1 * " 1 

fW 4 = 2 y cos a n y Sin 2 R m = -(AT + C 2a )S 2R 

n—1 m — 1 

Specific evaluation for design D 0 : 

P X QA ~ 9 

10. r_4_ 4 

General evaluation: 

K 

r xAA = = Slwl y <t>l M = Slw 2 xPxAA 

k= 1 


(195) 


(196) 


(197) 


(198) 


(199) 


(200) 


(201) 


( 202 ) 
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Using equations (58) gives 


K K 

Px AA = 4*1 Ai = Y1 '' OS ' C ° S2 ^ 

1 fc=l 

Specific evaluation for design D: 

If M 

Px AA = cos 2 a„ cos" /i',„ = — (TV + C 2a ){M + C 2 r) 

71=1 777=1 

Specific evaluation for design D 0 : 


P*AA 


L M (N + C. 2a ) 


(203) 


(204) 


(205) 


Evaluation of R, Matrix for y- and z- Axis Sensors 


By using the values of equations (58), the simplified elements of gradient vector <j> yc for the 
y-axis sensor are given by 


/» = 1 

f !/s = -cos a sin R = <f> XQ 

> 

(f> ya = —cos a cos R = <j> x 4 

4>y A = sin a = f XR 

and the elements of gradient vector <£, o for the z-axis sensor are given by 


(206) 


= 1 

/„ = —cos a cos R = <j> x 

<f>z Q = sin a = f xg 

4>z a = - cos a sin R = A 


(207) 


Equation sets (175), (206), and (207) show that the elements of vectors <f> yc and <j> % are 
permutations of vector <jj Xc . Since matrices P v and P;, are obtained from vectors <f> yc and <f> Zc) 
their rows and columns are permutations of matrix P t) and are the same permutations as those of 
vectors <f> yc and cf> Zc , respectively, relative to vector <j> Xc . Therefore, it follows from equation (127) 
that quadratic forms (jR r = qr^ = qR. and thus variance functions <j> Vx (z) = cr v ( z) = a Vz (z). 


Rj/ Matrix for y-Axis Sensor 

The elements of gradient matrix are obtained from equations (206) and (176) to (205) 
in terms of R a , as follows: 

r m = r *66 
r .%9 = P x bQ 
Pm = Pha 
Pm = r *6.9 - 



Similarly 


r V.9.9 — P X QQ 

P«sn — Pxqa 

Pvsa — Pxsa 


(209) 
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and 

Pm f! = P X AA 

PmA = P X SA > 
PvAA = r ®.9.9 ' 


R z Matrix for z-Axis Sensor 

The elements of gradient matrix R., are obtained from equations 
terms of as follows: 

r m — r m 

r zbS ~ P x bA ^ 

P*bn ~~ 1 x bs 
P%a = P*bn ■ 

Similarly 

'•,9.9 = P x AA ' 

Pz. 9SJ — Px. 9.4 * 

Pm A = P*C1A 

and 

Pma ~ r ®.9.9 
Pm . a z P x .mi > 

P Z AA — P-Hm ' 


(210) 

(207) and (176) to (205) in 
( 211 ) 

( 212 ) 

(213) 
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Appendix E 

Evaluation of Figure of Merit of Experimental Design 

The design figure of merit V for experimental design D is given by equation (53). The 
numerator of equation (53) contains integrals of cross products of the elements of gradient 
vector <f> c , which are now evaluated by using the parameter values of equations (58). 

The design figure of merit for the x-axis sensor is obtained from equations (21) to (26) 
and (58) as follows: 

I- = 1 


pA max /‘-ft i 




-‘nun -“mm 


/"I’m 
J ft™ in 


dR da = Ad A R 


(214) 


2 - fzjzs = sin a 


(‘Amax /‘ft] 


I x bS — 


-‘mm v -“nun 


/‘ftmax 

I sin 

d R-n-nr, 


adRda = — Ait A cos a 


(215) 


3. fz h 4>zn = -cos a sin it 


(‘A max /‘ftmaa 




/ 


cos a sin RdRda = A sin a A cos it 


(216) 


4- f n (t> XA = -cos a cos it 


rev 111 ax C-R Dai 




cos ce cos i? dR da = —A si n a A si n i? 


Amin J ftn: 


(217) 


5. fx s = sin 2 a 


max fftma 


^*.9.9 — 


/ 

^ft. 


sin a 


2 dit da = — A it I Aa — — A sin 2a 1 dR da 


(218) 


6. fz s <l>z Q — sin 2a sin it 


(‘Amax /‘ftma.x 


•**.9n — 2 


[ -*nun I'mm 


sin 2a sin RdRda = — — A cos 2a A cos it. (219) 


7. /*,</>*, = - —sin 2a cos it. 


^*.9A — 


^ f* °max f‘ ftn: 

2 J a...:.. J R.„. 


sin 2a cos RdRda = —A cos 2a A sin it 


( 220 ) 


‘min -“mm 


8. 4> XQ 4> xn = cos 2 a sin 2 it 


I XQS1 — 


(‘Am ax /‘ftmax 


' Anun ^ ft min 


1 


cos a sm 


1 


RdRda = — Aa + — A sin 2a ( A it — — A sin 2 it 


1 


( 221 ) 
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9. (f> XQ 4 ) x i — “ c °s 2 a sin 2 R 


'*nia.x r 1 ax 


IxaA 2 


cos 2 a sin 2 R dRda = — ( Act + — A sin 2a ) (A cos 2 R) (222) 

8 v 2 


10. (frxt&XA = cos 2 a cos 2 R 


-*max 1 1 l m ax 


^ x AA ~ 


1 


cos 2 a cos 2 RdR da = — I Aa + — A sin 2a 


1 


AT? + — A sin 2 R 


where 


Aa r r niax a , 


A R= 77, 


max 1 tmm 


Rn 


and 


A sin a = sin a max - sin a min 1 
A sin R = sin 7? max — sin 77,™,, J 

A sin 2a = sin 2a max — sin 2a min | 

A sin 2 R = sin 277 max — sin 2 77,™,, I 

Similar definitions apply for A cos a, A cos R , A cos 2a, and A cos 2 77. 
Define the following matrix where subscript x is omitted: 


4 = 


7 m 

hs 

IbQ 

4.4 


4.9 

7sq 

7.9.4 

4« 


7«n 

7(2.4 

I 

!*»■ 

7.9.4 

7(14 

7.4.4 _ 


It follows that 


(' rOc m a.x r Ri 11 ax 

V N = I q R { z) dx = / / 


4 4 


(223) 

(224) 

(225) 

(226) 


(227) 


-1 J.T 




c?77 da EE 1 ’ • 


(228) 


f=l i=i 

The figure-of-merit expression follows from equations (213), (214), and (228) as 

MNV n 


V = 


*66 


(229) 


The final expression applies to x - , y-, and 2 -axis sensors. 
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(<V CT *)’ V N' / 2 (a y Ja x ), V 


3.0 




(b) Calibrated from —45° to 45°. 

Figure 1. Normalized standard deviation of predicted output of single-axis AO A sensor without 
roll. 
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(MN) 1 / 2 (g Vx Ig x ), V (MN) ]/ 2 (a Vx /a x ), V 



(a) Calibrated from —30° to 30°. 



(b) Calibrated from —45° to 45°. 

Figure 2. Normalized standard deviation of predicted output of single-axis AOA sensor with roll. 
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Figure 3. Normalized standard deviation of predicted output of single-axis AOA sensor with roll 
for calibration points unequally spaced from — 30 c to 30°. 
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Figure 4. Normalized standard deviation of predicted output of single-axis AOA sensor with roll 
for calibration repeated at end points (±30°) and once at 0°. 



Figure 5. Normalized standard deviation of predicted output of single-axis AOA sensor with roll. 
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Pitch angle, deg 

Figure 6. Normalized standard deviation of inferred pitch angle of single-axis AO A sensor without 
roll for <f> a = 0°. 
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Pitch angle, deg 

(b) Q = 1°; A x = 20°; A y = 90°. 

Figure 9. Singularity loci of Jacobian matrix F, of x-y axis AOA eensor. 
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Pitch angle, deg 


(d) Q = 1°; A x = 90°; A y 


Figure 9. Concluded 







a lA a \\, a o/ G 



Roll angle, deg 

(a) cr x = 1; Q x = 1°; A x = 90°; a y = 1 ; fi v = 1°; A y = 0°. 


Figure 10. Normalized standard deviations of inferred pitch and roll angles of x-y axis AOA 
sensor. 
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Roll angle, deg 

(d) a x = 1; Q x = 1°; A x = 0°; a y = 10; Q :V = 0.1°; A y = 0°. 
Figure 10. Concluded. 
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(d)fi=l°;A, = 0 o ;A, = 0 o . 


Figure 11. Concluded 
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Pitch angle, deg 

Figure 13. Normalized standard deviations of inferred pitch and roll angles versus pitch angle for 
three-axis AOA sensor for cr y = a z = lOec,., Q x = Q y = fi x = 0.1°, and A x = 90°, A y = A z = 0°. 
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Figure 16. Normalized standard deviations of inferred pitch and roll angles versus roll angle for 
three-axis AOA sensor for a y = a, = lOer,., Q.*, = O i; = = 1°, and A x = 90°, A y = A z = 0°. 
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Figure 18. Normalized standard deviations of inferred pitch and roll angles versus roll angle for 
three-axis AOA sensor for = a z = <7* = 1, Q* = Q v = Q* = 1°, and A x = 90°, A y = A z = 0°. 




72 






gr/g v v CT r/7 c ' 




Figure 19. Normalized standard deviations of inferred pitch and roll angles versus pitch angle for 
three-axis AOA sensor for a y = a, = lOcr^., = Q ;/ = = 5°, and = 90°, A v = A z = 0°. 
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CTfl/CTy 



Figure 20. Normalized standard deviations of inferred pitch and roll angles versus roll angle for 
three-axis AOA sensor for a y = cr, = lOec,,, = O i; = Vl z = 5°, and A x = 90°, A y = A z = 0°. 
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Figure 22. Residuals of predicted output of single-axis AOA sensor without roll for six replications 
from —36° to 36°. 
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Figure 23. Errors of inferred pitch angles of single-axis AOA sensor without roll for six 
replications from —36° to 36°. 
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Figure 24. Residuals of predicted output of single-axis AOA sensor without roll for single-axis 
AOA sensor for six replications from —180° to 180°. 
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Figure 25. Errors of inferred pitch angles of single-axis AOA sensor without roll for six 
replications from —180° to 180°. 
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Figure 26. Residuals of predicted output of single-axis AOA sensor without roll for six replications 
and four-point tumble test. 
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Pitch angle, deg 

Figure 27. Errors of inferred pitch angle of single-axis AOA sensor without roll for six replications 
and four-point tumble test. 
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(a) Without temperature correction. 


Figure 28. Residuals of predicted output of single-axis AOA sensor with roll for six replications 
from —30° to 30°. 
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(b) With temperature correction. 
Figure 28. Concluded. 
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Figure 29. Errors of inferred pitch angle of single-axis AOA sensor with roll for six replications 
from —30° to 30°. With temperature correction. 
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Figure 30. Errors of inferred pitch angle of single-axis AOA sensor with roll for one replication 
from —30° to 30°. With temperature correction. 
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Roll angle, deg 

Figure 31. Residuals of predicted output versus roll angle of single-axis AOA sensor with roll for 
six replications from —180° to 180°. With temperature correction. 
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Figure 33. Errors of inferred pitch angle versus roll angle of single-axis AOA sensor with roll for 
one replication from —180° to 180°. With temperature correction. 


88 



Predicted output error, deg 



Pitch angle, deg 


Figure 34. Residuals of predicted output of single-axis AOA sensor 2 for six replications from 
—30° to 30°. With temperature correction. 
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Figure 35. Residuals of predic ted output of single-axis AOA sensor with roll for fractional design 
and six replications from —30° to 30°. With temperature correction. 
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Figure 36. Residuals of predicted output of single-axis AOA sensor with roll that were 
recomputed by using parameters estimated from fractional design. With temperature 
correction. 
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Figure 38. Errors of inferred pitch angle of single-axis AOA sensor with roll for four replications 
from —180° to 180°. With temperature correction. 
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Figure 39. Residuals of predicted output of singl^axis AOA sensor with roll for one replication 
from —180° to 180°. With temperature correction. 
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Figure 40. Residuals of predicted output versus roll angle of single-axis AOA sensor with roll for 
four replications from —180° to 180°. With temperature correction. 
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Figure 41. Residuals of predicted output versus roll angle of single-axis AOA sensor with roll for 
one replication from —180° to 180°. With temperature correction. 
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Figure 42. Residuals of predicted output of single-axis AOA sensor 2 with roll for six replications 
from —180° to 180°. With temperature correction. 
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Figure 44. Con 
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Figure 45. Errors of inferred pitch angles of three-axis AOA package with roll for one replication 
from —90° to 90°. With temperature correction. 
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Figure 46. Errors of inferred roll angles of three-axis AOA package with roll for one replication 
from —90° to 90°. With temperature correction. 
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Figure 47. Predicted output residuals of three-axis AOA package with roll for six replications 
from —180° to 180°. With temperature correction. 
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(b) y-axis sensor. 
Figure 47. Continued. 
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Figure 47. Concluded. 
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Figure 48. Errors of inferred pitch angles of three-axis AOA package with roll for one replication 
from —180° to 180°. With temperature correction. 
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Figure 50. Errors of predicted output residuals of x-, y-, and 2 -axis sensors of three-axis 
AOA package with roll for four-point tumble test with six replications. With temperature 
correction. 
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Figure 51. Errors of inferred pitch and roll angles of three-axis AOA package with roll for six- 
point tumble test with six replications. With temperature correction. 
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Figure 52. Predicted output residuals of three-axis AOA package with roll calculated by using 
parameters estimated from six-point tumble test. With temperature correction. 
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(b) y - axis sensor. 
Figure 52. Continued. 
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Figure 53. Predicted output residuals of .c-axis sensor of three-axis AOA package with roll for 
fractional design with six replications. 
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Figure 54. Predicted output residuals of a;-axis sensor of three-axis AOA package with roll 
calculated by using parameters estimated from fractional design. 
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